Draft version February 28, 2013 

Preprint typeset using I^T^X style cmulateapj v. 5/2/11 



A KINEMATIC STUDY OF THE ANDROMEDA DWARF SPHEROIDAL SYSTEM 

Michelle L. M. Collins 1 ' 2 , Scott C. Chapman 2,3 , R. Michael Rich 4 , Rodrigo A. Ibata 5 , Nicolas F. Martin 5 ' 1 , 
Michael J. Irwin 2 , Nicholas F. Bate 6 , Geraint F. Lewis 6 , Jorge Penarrubia 7,8 , Nobuo Arimoto" 10 , Caitlin M. 
Casey 11 , Annette M. N. Ferguson 8 , Andreas Koch 12 , Alan W. McConnachie 13 , Nial Tanvir 14 

Draft version February 28, 2013 

ABSTRACT 

We present a homogeneous kinematic analysis of red giant branch stars within 18 of the 28 An- 
dromeda dwarf spheroidal (dSph) galaxies, obtained using the Keck I LRIS and Keck II DEIMOS 
spectrographs. Based on their g — i colors (taken with the CFHT MegaCam imager), physical posi- 
tions on the sky, and radial velocities, we assign probabilities of dSph membership to each observed 
star. Using this information, the velocity dispersions, central masses and central densities of the dark 
matter halos are calculated for these objects, and compared with the properties of the Milky Way 
dSph population. We also measure the average metallicity ([Fe/H]) from the co-added spectra of 
member stars for each M31 dSph and find that they are consistent with the trend of decreasing [Fe/H] 
with luminosity observed in the Milky Way population. We find that three of our studied M31 dSphs 
appear as significant outliers in terms of their central velocity dispersion, And XIX, XXI and XXV, 
all of which have large half-light radii (> 700pc) and low velocity dispersions (a v < 5kms _1 ). In 

addition, And XXV has a mass-to-light ratio within its half-light radius of just [M/L]haif = 10.3^g'!j!, 
making it consistent with a simple stellar system with no appreciable dark matter component within 
its la uncertainties. We suggest that the structure of the dark matter halos of these outliers have 
been significantly altered by tides. 

Subject headings: dark matter — galaxies: dwarf — galaxies: fundamental parameters — galaxies: 
kinematics and dynamics — Local Group 



1. INTRODUCTION 

The underlying nature of the dark matter halos of 
dwarf spheroidal galaxies (dSphs) has garnered signif- 
icant attention from the scientific community over the 
past decade. The goal of recent observational studies of 
these objects has been to make critical tests of structure 
formation scenarios, particularly focusing on the viability 
of the canonical ACDM model. There is the long stand- 
ing issue of the relative dearth of these faintest of galax- 
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ies observed surrounding nearby galaxies when compared 
with the number of dark matter subhalos produced in 
N— body simulations, which is referred to as the "miss- 
ing s atellite" problem (jKlvpin et all 119991 : iMoore et all 
1999). The extent to which this mismatch is considered 
problematic has decreased over recent years as both theo- 
rists and observers have sought to reconcile the simulated 
and observable Universe. From a modelling point of view, 
one does not expect stars to be able to form within all 
dark matter subhalos seen in simulati ons, and at a cer- 
tain mass limit (Vmax 5, lSk ms" 1 , see iPenarrubia et "all 
l2008al IKoposov et al.l |2009) , star formation is unable 
to proceed. Thus, there is a lower limit placed on 
galaxy formation. This mass limit is also tied to feed- 
back processes that can remove th e baryonic reservoirs 
required for star formation (e.g. , IBullock et al.l 
Somervilld l200l iKTavtsovl [2010: IBullock et all 



2000; 



2010: 



Nickerson et al.l l2011t iKazantzidis et al.l 120111 1 This 
would imply that only the most massive subhalos seen 
in simulations are able to form and retain luminous pop- 
ulations. Observers have also attempted to quantify cur- 
rent survey completeness and radial selection effects to 
account for the num ber of satellites we are not currently 
able t o detect (e.g . . IKoposov et al.l 12008: Tol lerud et al.l 
120081: iWalsh et all I2009D . These studies suggest that 
there are of order a few hundred satellites within the 
Milky Way's virial radius that we have yet to detect. 

The high dark matter dominance of dSph galaxies 
also singles them out as objects of interest. With to- 
tal dynamic mass-to-light ratios of [M/L] ~ 10 — 1000s 
and half-light radii of r^aif ~ 10 — 1000 pc, they are 
ideal systems with which to probe the inner density pro- 
files of dark matter halos. Recent imaging and spec- 
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troscopic observations of these objects within the Lo- 
cal Group have shown that, despite spanning approxi- 
mately 5 decades in luminosity, the dSphs of the MW 
shar e a common mass scale and a universal density pro- 
file (I Strigari et alJl2?M [Walker et alJl2009ht IWolfeTai] 
12010( 1. With a kinematic resolution of a few 10s of par- 
sees, these objects allow us to start addressing the ques- 
tion of whether the central regions of these halos follow 
cuspy de nsity profiles as pred icted by cosmological sim- 
ulations (|Navarro et al.lfl997t ). or constant density cores, 
similar t o what is observed in low surface brightness 
galax i es dBlais-Ouellette et al.l 120011: Ide Blok fe Bosnia] 
200 1 IdeBlok et all 120031: Ide Blokl 120051: iSwaters et al.l 
200a iKassin et alj|2006t iSpano et al.ll2008[ ). From stud- 
ies of brighter dSphs , such as Sculptor and Fornax 
dWalker fc Penarrubial 120111 : lAmorisco fe Evansl 120121 : 
Uardel fcG cbhardtl l2012h it appears that their halo den- 
sity profiles are also inconsistent with hosting central 
cusps. It is possible that these objects originally formed 
with cuspy density profiles, and that these have been sub- 
sequently modified by baryonic feedback. If this is truly 
the case, one should be able to observe cuspy profiles in 
the fainter dSph population (|Zolotov et al.ll2012i ) as these 
do not contain enough baryons to drive this change in the 
dark matter density profile. Perhaps the only way to gain 
further insight into this contentious issue is by measuring 
the kinematics for large numbers of stellar tracers within 
these objects and analysing them with detailed models 
that do not make ab initio assumptions about the un- 
derlying density profiles, or the velocity anisotropy of 
bo th the dark matter and stars , suc h as those employed 
bvlWalker fc Peharrubial (|201lD and Uardel fc Gebhardti 
(|2012D . 

To date, the majority of studies involving the de- 
tailed kinematics of dSphs have revolved largely around 
those belonging to the Milky Way, as these are nearby 
enough that we can measure the velocities of their 
member stars to a high degree of accuracy. How- 
ever, there are currently only ~ 25 known MW dSphs, 
with luminosities ranging from 10 2 - 1O 7 L . For the 
very faintest, some controversy remains as to whether 
they are massively dark-matter dominated ( see e.g., 
iNiederste-Ostholt et al.ll2009t ISimon et al.ll201lb . but al- 
most all of them have been sh own to be consistent 
with the universal mass profiles of iWalker et~aH (j2009al) 
and IWolf et all (|2010l ) . One notable exc eption to this 
is the Hercules object ([Aden et al.l 12009( 1 , which some 
have arg ued is currently unde rgoing significant tidal dis- 
ruption (|Martin fc Jinl 120101 ) . Andromeda represents 
the only other system for which comparable kinematic 
analyses can be performed. M31 now has 28 dSph 
companions known, whose luminosities range from ~ 
10 4 — 10 8 L Q , the majority of which have been dis- 
covered by the C FHT Pan-Androm e da Archaeological 
Survey (PAndAS iMartin et al.l 120061: llbata et "all 120071: 
Irwin et al.ll2008HMcCoimachie et"aHl2008l; iMartin et al.l 
20091: iRichardson et al.l 1201 ID . The relatively brighter 
lower bound for the luminosities of M31 dSphs com- 
pared to the MW is a detection limit issue, rather than 
a sign of differing stellar populations (Martin et al. 
2013, in pr ep). It has been noted by a number of au- 
thors fe.g..lMcConnachie fc Irwinll2006bl: iTollerud et al.l 
l20ll IMcConnachid 120121) that for the brighter dSphs 
(My < —8), those belonging to M31 are 2-3 times more 



extended in terms of their half-light (rhaif) and tidal (rt) 
radii compared with the MW. In these papers, the under- 
lying cause of this discrepancy was not identified, but it 
has been argued that it could be an effect of environment, 
with the mass distribution of th e host playing an impor - 
tant role. Subsequent work bv Br asseur et ahl ()2011a( ). 
who included the fainter, non-classic M31 dSphs for the 
first time, showed that statistically, the relationship be- 
tween size and luminosity for dSphs in the MW and An- 
dromeda are actually largely consistent with one another, 
however there remain a number of significantly extended 
outliers within the Andromcdcan system (e.g., And II, 
And XIX, rhaif ~ 1.2 kpc and 1.5 kpc respectively), and 
the scatter in this relationship is large (up to an orde r of 
magnitude at a given luminosity. lMcConnachie]|2012f) . 

Wo rking from the iMcConnachie fc Irwinl (|2006b| ) re- 
sults, iPenarrubia et al.l (|2008aj ) modelled the expected 
velocity dispersions for the M31 dSphs, assuming that 
all dSph galaxies are embedded within similar mass dark 
matter halos. A robust prediction of their modeling was 
that, given the larger radial extents, the dSphs of M31 
should be kinematically hotter than their MW counter- 
parts by a factor of ~ 2. At the time of writing, they 
had only 2 measured velocity dispers ions for the M31 
dSphs, those of And II and And IX (|C6te et all 119991 : 
IChapman et al.ll2005f). New studies of the kinematics of 
M31 dSphs (iCollins et aLll2010l l20Tll iKalirai et al.l l2010l: 
ITollerud et al.ll2012tlChapman et alil2012D have dramat- 
ically increased the number of systems with a measured 
velocity dispersion, and have shown that instead of be- 
ing kinematically hotter, these systems are either very 
similar to, or in a number of cases (e.g., And II, And 

XII, And XIV, And XV and And XXII), colder than 
their MW counterparts. In particular, a significant re- 
cent kinematic study of 15 M31 dSph companions us- 
ing the Keck II DEIMOS spectrograph by the Spectro- 
scopic and Photom etric Landscape of An dromeda's Stel- 
lar Halo (SPLASH. ITollerud et al.ll20Tl concluded that 
the M31 dSph system largely obeys very similar mass- 
size-luminosity scalings as those of the MW. However, 
they also identified 3 outliers (And XIV, XV and XVI) 
that appear to possess much lower velocity dispersions, 
and hence maximum circular velocities, than would be 
expected for these systems. Such a result suggests that 
there are significant differences in the formation and/or 
evolution of the M31 and MW dSph systems. 

To investigate this further, our group has been sys- 
tematically surveying the known dSphs of M31 with the 
Keck I LRIS and Keck II DEIMOS spectrographs, and 
have obtained kinematic data for 18 of the 28 galaxies. 
In this paper, we present new spectroscopic analysis for 
the 11 dSphs, Andromeda (And) XVII, And XVIII, And 
XIX, And XX, And XXI, And XXIII, And XXIV, And 
XXV, And XXVI, the tidally disrupting And XXVII, 
And XXVIII and And XXX (Cassiopeia II) using an al- 
gorithm we have developed that implements a probabilis- 
tic method of determining membership for each galaxy. 
In addition we re-analyze the kinematics of 6 dSphs that 
our group has previously observed (And V, VI, XI, XII, 

XIII, XXII) using this method with the aim of provid- 
ing a homogeneous analysis of all dSphs observed by our 
group to date. We also provide the individual stellar 
velocities and properties for every star observed in our 
dSphs survey, allowing us to present a large catalog of 
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stellar kinematics that will be of interest to those study- 
ing dSph systems and Milky Way-like galaxies, whether 
observationally or theoretically. 

The outline of this paper is as follows. In § 2 we discuss 
the relevant observations, data reduction techniques. In 
§ 3 we outline our algorithm for the classification of mem- 
ber stars. In § 4 we present an analysis of our new kine- 
matic datasets. In § 5 we report on the masses and den- 
sities of the dark matter halos of the M31 dSphs, com- 
paring them to those of the MW dSphs. In §6 we report 
on the metallicities, [Fe/H], of the M31 dSphs as mea- 
sured from the co-added spectra of their member stars. 
Finally, we conclude in § 7. 

2. OBSERVATIONS 
2.1. Photometry and target selection 

The PAndAS survey (jMcConnachie et al J [2009( 1. con- 
ducted using the 3.6 metre Canada France Hawaii Tele- 
scope (CFHT), maps out the stellar density of the disc 
and halo regions of the M31-M33 system over a pro- 
jected area of ~ 350deg 2 (~ 55,000kpc 2 ), resolving in- 
dividual stars to depths of g = 26.5 and i = 25.5 with a 
signal to noise ratio of 10, making this survey the deep- 
est, highest resolution, contiguous map of the majority of 
the extended stellar halo of an L* galaxy to date. Each 
of the 411 fields in this survey (0.96 x 0.94 deg 2 ) has 
been observed for at least 1350s in both MegaCam g 
and i filters, in < 0.8" seeing. This survey was ini- 
tiated following two precursor surveys of the M31 sys- 
tem, the first of which surveyed the central ~ 40 deg 2 
conducted with the 2.5 metre Isaac N ewton Telescope 
(jFerguson et all [20021: H rwin et ail 120051 ). and revealed a 
wealth of substructure in the Andro meda stellar halo , 
including the giant southern stream (jlbata et al.l 120011 ) . 
To better understand this feature, and to probe deeper 
into the M31 (and M33) stellar halo, a survey of the 
south west qu adrant of the M31 halo was initiated us- 
ing the CFHT (jlbata et al.ll2007f) . and revealed yet more 
substructure, includin g the arc like stream Cp and Cr 
([Chapman et a l. 2008) and a number of dwarf spheroidal 
satellites ([Martin et all 12006( 1. This CFHT survey was 
then extended into the full PAndAS project. For de- 
tails of the proce s sing a nd reduction of these data, see 
iRichardson eTaTl (|20lTf ). This survey has introduced 
us to a wealth of stellar substructure, debris and glob- 
ular clusters within the Andromeda- Triangulum sys- 
tem. In addition, it has led to the discovery of 17 
dSphs. These objects were detected in the PAndAS sur- 
vey maps as over-densities in matched-filter surface den- 
sity maps of metal po or red giant branch ( RGB) stars 
and we re presented in IMartin et al.l (120061): llbata et aLl 
(120071): llrwin et all (120081): iMcConnachie et al.1 (12 008): 
IMartin et al.1 ([2009( 1 and IRichardson et al.l ([20111 ). We 
briefly summarise the photometric properties of all dSphs 
discussed within this paper in Table [TJ 

For the majority of these objects, the PAndAS 
dataset formed the basis for our spectroscopic target 
sel ection. Using the color s election boxe s pres ented 
in IMcConnachie et al.l (120081) : IMartin et ahl ([2009D and 
IRichardson et al.l ~ 201lD . we isolated the RGBs of each 
dSph, then prioritised each star on this sequence depend- 
ing on their color, i-band magnitude, and distance from 
the centre of the dSph. Stars lying directly on the RGB, 



with 20.3 < io < 22.5 and distance, d < 4rh a if (where 
'"half is the half-light radius, measured on the semi-major 
axis of the dSph) were highly prioritised (priority A), 
followed by stars on the RGB within the same distance 
from the centre with 22.5 < io < 23.5 (priority B). The 
remainder of the mask was then filled with stars in the 
field with 20.3 < i < 23.5 and 0.5 < g - i < 4 (priority 
C). In general, it is the brighter, higher priority A stars 
that ultimately show the highest probability of member- 
ship. We also use the PAndAS photometry to help us 
determine membership of the dSph (discussed in §[3]). 

And XXVIII is the one dSph in our sample that is 
not covered by the PAndAS survey. For this object, tar- 
get selection followed an identical methodology tot hat 
detailed above, but used photometry from the 8th data 
release of the Sloan Digital Sky Survey (SDSS-III). De- 
tails of the obse rvations and analys is of this photometry 
can be found in lSlater et al.l ([201 ID . 

2.2. Keck Spectroscopic Observations 

The DEepTmaging Multi-Object Spectrograph 
(DEIMOS), situated on the Nasmyth focus of the Keck 
II telescope is an ideal instrument for obtaining medium 

o 

resolution (R~ 1.4 A) spectra of multiple, faint stellar 
targets in the M31 dSphs. The data for the dSphs within 
this work were taken between Sept 2004 and Sept 2012 
in photometric conditions, with typical seeing between 
0.5 — 1". Our chosen instrumental setting covered the 

o 

observed wavelength range from 5600-9800 A and em- 
ployed exposure times of 3x20 minute integrations. The 
majority of observations employed the 1200 line/mm 
grating, although for 4 dSphs (And XI, XII, XIII 
and XXIV) the lower resolution 600 line/mm grating 

(R ~ 3.8 A FWHM) was used. The spectra from both 
setups typically possess signal-to-noise (S:N) ratios of 

o 

> 3 A -1 for our bright targets (i < 22.0), though some 
of our fainter targets fall below this level. Information 
regarding the spectroscopic setup and observations for 
each dSph are displayed in Table [2] 

The resulting science spectra are r educed usin g a cus - 
tom built pipeline, as described in llbata et al.l ([2011( 1. 
Briefly, the pipeline identifies and removes cosmic rays, 
corrects for scattered light, performs flat-fielding to cor- 
rect for pixel-to-pixel variations, corrects for illumina- 
tion, slit function and fringing, wavelength calibrates 
each pixel using arc-lamp exposures, performs a 2- 
dimensional sky subtraction, and finally extracts each 
spectra - without resampling - in a small spatial region 
around the target. This results in a large set of pixels for 
each target, each of which carries a flux and wavelength 
(with associated uncertainties) plus the value of the tar- 
get spatial profile at that pixel. We then derive velocities 
for all our stars with a Bayesian approach, using the Ca II 
Triplet absorption feature. Located at rest wavelengths 

o 

of 8498, 8542 and 8662 A, these strong features are 
ideal for determining the velocities of our observed stars. 
We determine the velocities by using an Markov Chain 
Monte Carlo procedure where a template Ca II spectrum 
was cross-correlated with the non-resampled data, gen- 
erating a most-likely velocity for each star, and a likely 
uncertainty based on the posterior distribution that in- 
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TABLE 1 

Details of the structura l properties for each dwarf, as derived from 
CFHT photometry byIIbata et AL.I (120071); ILetARteet AL.I (T2009I); 
Irwin et al. (2008); McConnachie et al. (2008); Martin et al. (2009); 
Collins et al. (2010. 2011); Brasseur et al. (2011a); Richardson et al. 
(120111) and [McConnachTeI l|2012l) . updated for distances presented in 
IConn et al.I (120121 ). The p hotometry for And VI is derived from Subaru 
SuprimeCam photometry (ICollins et al.II2011|), and the values f or And 
XXVIII come from SDSS photometry IISlAter et al.1120111) . 
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corporates all the uncertainties for each pixel. Typically 
our velocity uncertainties lie in the range of 5 — 15 km s . 
Finally, we also correct these velocities to the heliocentric 
frame. 

2.2.1. Telluric velocity corrections 

With slit-spectroscopy, systematic velocity errors can 
be introduced if stars are not well aligned within the 
centre of their slits. Such misalignments can result from 
astrometric uncertainties or a slight offset in the position 
angle of the mask on the sky. For our astrometry, we take 
the positions of stars from PAndAS photometry, which 
have an int ernal accuracy of ~ 0.1" and a global accuracy 
of ~ 0."25 (|Segall et alJl2007h . This can translate to ve- 
locity uncertainties of up to ~ 15kms _1 for ourDEIMOS 
setup. In previous studies, authors have tried to correct 
for this effect by cross correlating their observed sp ectra 
with telluric absorption features (e. g.. ISohn et alj|2007t 
Simon fc Gehal 120071: iKalirai et al.l 120101 : ICollins et al l 
20101: iTollerud et al.ll2012l ). These atmospheric absorp- 
tion lines are superimposed on each science spectrum, 
and should always be observed at their rest-frame wave- 
lengths. Thus, if one is able to determine the offset of 
these features, shifts caused by misalignment of the sci- 
ence star within the slit can be corrected for, and this 
can be applied on a slit by slit basis. The strongest of 
these features is the Fraunhofer A-band, located between 

o 

7595-7630 A. An example of this feature is shown in the 
left panel of Fig. [T] 

While we believe this correction is robust in the high 
S:N regime, we argue against applying this correction 
in studies of Andromedean satellites, where S:N is often 



quite low (typically less than 8-10 A -1 ) for 1 hour ob- 
servations of faint (i > 21) RGB stars. We find that 
when we compute this offset for all stars within our sam- 
ple, those with high S:N tend to cluster within a few 
kms -1 of the average telluric correction found for the 

o 

mask. As the S:N decreases below about 10 A -1 , the 
scatter about this mean value increases dramatically, as 
do the uncertainties computed for each individual cor- 
rection. This is because the telluric feature is a single, 
very broad and asymmetric feature. It is therefore easy 
in the noisy regime for the cross correlation routine to 
misalign the template and science spectrum whilst still 
producing a high confidence cross-correlation maximum. 
We show this effect explicitly in the right hand panel of 
Fig. [T] where we plot the deviation of the telluric correc- 
tion for every star within our sample from the average 
correction determined for the spectroscopic mask it was 
observed with (vtei— < Vtei >) as a function of S:N. The 
cyan points represent the individual data, and the large 
black points represent the median of all points within 

o 

1 A bins in S:N. The error bars represent the dispersion 
within each bin. It is plainly seen that the median value 
for each bin is consistent with zero (i.e., the median), and 
that the dispersion increases with decreasing S:N. If we 
were to apply these velocity corrections to all our stars, 
it is probable that we would merely increase the velocity 
uncertainties rather than reducing them. 

For this reason, we take a different approach. Using 
solely the telluric velocity corrections of stars from each 
observed mask whose spectra have S:N> 7, we measure 
(a) the average telluric correction for the mask and (b) 
the evolution of the telluric correction as a function of 
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TABLE 2 

Details of spectroscopic observations 
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C«0,J2000 


<5(),J2000 


P.A. 


Exp. (s) 


Nxarg 


NMem 


And V 


Ifi A no- 200Q 


LRIS 


831 /8200 


o 

3.0 A 


01 -1 0-1 8 21 


+47-37-53 3 


0° 


3600 


50 


15 


And VI 


1 7-1 9 Sent 2009 


DETMOS 


1200 


o 

1.3 A 


23-51 -51 4Q 


4-24-34-57 


0° 


5400 


113 


43 


And XI 


93 Sonf 900fi 


DFTMOS 

±y nil ivi vj o 


600 


o 

3.5 A 


UU.'iU.^O.UO 


TOO . 1U .ZiO.O 


0° 


3600 


33 


5 


And XTT 


2 1 — 2 S Sent 2006 


DETMOS 

i /ill ivi v_y o 


600 


o 

3.5 A 


00-47-39 SQ 


4- 34-99-98 6 


0° 


3600 


49 




And XTTT 


23 Sent 200fi 


DETMOS 

1 / 1 t 1 IVI V_y O 


600 


o 

3.5 A 


00-52-00 22 


4-32-59-16 2 


0° 


3600 


46 


4 


And XVTT 


2fi Sent 201 1 


DETMOS 

i / i 1 1 ivi v_y o 


1200 


o 

1.3 A 


00-37-51 09 


4-44-17-^1 9 


280° 


3600 


119 




And XVTTT 


Q Sent 2010 


DETMOS 

1 / 1 f 1 IVI V_7 O 


1200 


o 

1.3 A 


00-02-14 50 


4-45-05-20 


90° 


3600 


73 


4 


And XTXa 


2fi Sent 201 1 


DETMOS 

i / i 1 1 ivi v_y o 


600 


o 

3.5 A 


00-1 Q-45 04 


4-35-05-28 8 


270° 


3600 


107 


15 


And XTXh 


2fi Sent 201 1 


DETMOS 

1 / 1 f 1 IVI V_7 O 


1200 


o 

1.3 A 


00-1 Q-30 88 


4-35-07-34 1 


0° 


3600 


108 


9 


And XX 


q Sent 2010 


DEIMOS 


1200 


o 

1.3 A 


00-07-30 69 


+ 35-08-02 4 


90° 


3600 


85 


1 


And XXT 


2fi Sent 201 1 


DETMOS 

1 / 1 J 1 IVI v_y O 


1200 


o 

1.3 A 


23-54-47 70 


4-42-28-33 6 


180° 


3600 


157 


32 


And XXTTa 


23 Sent 2009 


DETMOS 

1 / 1 J 1 IVI v_y O 


1200 


o 

1.3 A 


01 -27-52 37 


4-28-05-22 3 


90° 


1200 


93 


4 


And XXTTh 


10 Sent 2010 


DETMOS 

1 / 1 J 1 IVI v_y O 


1200 


o 

1.3 A 


01 -27-52 37 


4-28-05-22 3 


0° 


3600 


73 


f; 


And XXTTTa 

iLllLi. Aj\lll(l 


q Q Pn t 201 


DETMOS 

1 / 1 J 1 IVI v_y O 


1200 


o 

1.3 A 


01 -29-1 8 1 8 


4- 38-43-50 4 

TOO . rtO . U*J . 


315° 


3600 


196 


21 


And XXIIIb 


10 Sept 2010 


DEIMOS 


1200 


o 

1.3 A 


01:29:21.87 


+38:44:58.7 


245° 


3600 


189 


18 


And XXIVa 


9 Sept 2010 


DEIMOS 


1200 


O 

1.3 A 


01:18:32.90 


+46:22:50.0 


30° 


3600 


192 


1 


And XXIVb 


31 May 2011 


DEIMOS 


600 


3.5 A 


01:18:32.90 


+46:22:50.0 


0° 


2700 


115 


11 


And XXV 


10 Sept 2010 


DEIMOS 


1200 


O 

1.3 A 


00:30:01.88 


+46:50:31.0 


90° 


3600 


183 


26 


And XXVI 


10 Sept 2010 


DEIMOS 


1200 


O 

1.3 A 


00:23:41.42 


+47:54:56.8 


90° 


3600 


179 


6 


And XXVII 


10 Sept 2010 


DEIMOS 


1200 


O 

1.3 A 


00:37:31.93 


+45:23:55.4 


45° 


3600 


131 


8 


And XXVIII 


21 Sept 2012 


DEIMOS 


1200 


O 

1.3 A 


22:32:32.71 


+31:13:13.2 


140° 


3600 


102 


17 


Cass II 


28 Sept 2011 


DEIMOS 


1200 


O 

1.3 A 


00:37:00.84 


+49:39:12.0 


270° 


3600 


156 


8 



Fraunhofer A Band 



3 
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.20 



10 
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-20 
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FlG. 1. — Left: Keck II DEIMOS spectrum of a bright foreground star, centered on the region of the strong telluric absorption feature, 
the Fraunhofer A band. Right: The deviation of the telluric correction for every observed star within our sample from the average 
correction determined for the spectroscopic mask (v te i — < vt e i >) as a function of S:N. Cyan points represent the individual data, and the 

o 

large black points represent the median of all points within 1 A bins in S:N. Error bars show the standard deviation within the bin. The 
median value for each bin is consistent with zero (i.e., the median), and the dispersion increases with decreasing S:N, arguing against using 
this correction in the low S:N regime. 
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mask position. In this way, we can track any gradient 
in our measurements that could be caused by e.g., ro- 
tational offsets in our mask. In general, we find these 
corrections to be slight. The average measured offset 
across all our masks is 3.8 km s" 1 (ranging from between 
— 3.4kms and +10.6kms _1 ). The measured gradients 
are very slight, resulting in an average end-to-end mask 
difference of 2.6 km s -1 , with a range of 0.1 — 7.2 km s -1 , 
typically within our measured velocity uncertainties. 

3. A PROBABILISTIC DETERMINATION OF MEMBERSHIP 

Determining membership for Andromedcan dwarf 
spheroidals is notoriously difficult in the best of cases. 
We only possess information about the velocity, CMD 
position, distance from the centre of the dSph and spec- 
troscopic metallicity (although this carries large uncer- 
tainties of > 0.3 dex for individual stars). Depending on 
the systemic velocity of the dSph, we must try to use 
these properties to distinguish the likely members from 



either Milky Way halo K-dwarfs (vh e 



> 



-150 kms 



or M31 halo giants (vhalo 



-300 kms \ a v halo 



90 kms 1 , Chapman ct al. 2006). In the case of Galactic 
contamination, our spectra also cover the region of the 

o 

Na I doublets (~ 8100 A). As this feature is dependent 
on the stellar surface gravity, it is typically stronger in 
dwarf stars than in giants. However, there is a significant 
overlap between the two, especially in the CMD color re- 
gion of interest for Andromedcan RGB stars. In the past, 
groups have focused on making hard cuts on likely mem- 
bers in an attempt to weed out likely contaminants (e.g., 
based on their distance from the centre of the galaxy or 
their velocity. iChapman et all 120051: iCollins et ail 120101 
120111: iKalirai et al .||2010D, bu t such 'by eye' t echni ques 
are not particularly robust. iTollerud et al.l (|2012l ) re- 
cently presented an analysis of a number of M31 dSphs 
where they used a more statistical method to ascertain 
likely membership, using the distance from the centre of, 
and position in the CMD of stars targeted within their 
DEIMOS masks. Here we employ a similar technique 
that will assess the probability of stars being members of 
a dSph based on (i) their position on the CMD; (ii) their 
distance from the centre of the dSph and, in addition; 
(hi) their position in velocity space, giving the likelihood 
of membership as: 



Pi OC P C MD X Pdist X P ve 



(1) 



Below, we fully outline our method, and implement a 
series of tests to check it can robustly recover the kine- 
matics of the M31 dSphs. 

3.1. Probability based on CMD position, Pcmd 

For the first term, Pcmd, we are interested in where a 
given star observed in our mask falls with respec t to the 
RGB of the observed dSph. ITollerud et ail ([20121 ) use the 
distance of a given star from a fiducial isochronc fit to 
the dwarf photometry to measure this probability. Here, 
we determine this value from the data itself, rather than 
using isochrones. Using the PAndAS CFHT photometry, 
we construct a normalised Hess diagram for the central 
region (i.e., within 2 x Hiaif) of the dSph, and one of 
a surrounding 'field' comparison region. By combining 
these two Hess diagrams, we can then map both the color 



distribution of the dSph and that of our contaminating 
populations. We use both directly as probability maps, 
where the densest region would have a value of 1. 

So we are not dominated by shot noise of sparsely pop- 
ulated regions of the CMD, we use only the region of 
the RGB. We do this by assigning a generous bounding 
box around the RGB as seen in the left hand panel of 
Fig. [2] where we display the PAndAS CMD for the well 
populated And XXI RGB. We have zoomed in on the 
region for which DEIMOS observations with reliable ve- 
locities can be obtained, e.g., i < 23.5. The bounding 
box is shown with red dashed lines. Anything that falls 
outside this region is therefore assigned a probability of 
Pcmd = 0. The resulting probability map for And XXI 
is shown in the right hand panel. Red points represent all 
DEIMOS stars that have Pcmd > 10~ 6 , while the blue 
stars show stars from the DEIMOS mask that are far re- 
moved from the And XXI RGB, and thus not considered 
to be members. 

3.2. Probability based on distance position, Pdist 

The second term in our probability function, Pdist 
can be easily determined from the known radial profile 
of the dSphs. The half-light ra dii of all these objects 



are known and can b e found in McConnachic & Irwin 
(|2006al): iZucker et all (T2004. 20 0~7ft: iMcConnachie et al 



2008); M artin et al.1 (120091 ): ICollins et al.l (l2010l I20T1 
Richardso n et al. l (|2011[ ). \¥e also know that their den- 
sity profiles are well represented by a Plummer profile 
with a scale radius of r p = rhaif • Therefore, we can de- 
fine th e probability fu nction as a normalised Plummer 
profile (|Plummerlll911l ). i.e.,: 



Pdist — 



7rr2[l + (r/r p ) 2 



(2) 



The above equation assumes that the systems we are 
studying are perfectly spherical. While the majority of 
these systems are not significantly elliptical, it is impor- 
tant to consider the effect of any observed deviations 
from sphericity. We therefore modify r p based on a given 
stars angular position with respect to the dwarfs major 
axis, 0i, such that: 



niaif(l - e) 
1 + ecosOi 



(3) 



wher e e is the measured e llipticity of the dSph as taken 
from lMcConnachid (|2012l ). 

3.3. Probability based on velocity, P ve i 

The final term, P ve i, contains information about the 
likelihood of a given star belonging to a kinematic sub- 
structure that is not well described by the velocity pro- 
files of either the MW halo dwarfs or the Andromeda 
halo giants, both of which are determined empirically 
from our DEIMOS database of > 20000 stars. From an 
analysis of these stars (selecting out any obvious sub- 
structures and the Andromedean disc) we find that the 
M31 halo is well approximated by a single Gaussian 
with a systemic velocity of v r ^haio — —308.8 kms -1 and 
o~v.haio = 96.3 kms -1 , giving a probability density func- 
tion for a given star with a velocity Vi and velocity un- 
certainty of v err ,i of: 
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Fig. 2. — Left: The PAndAS CMD for And XXI, showing all stars within 2 X r half of the dwarf centre. The red dashed lines show 
the bounding box we use to construct our probability map for the RGB of the dSph. Right: A background normalized Hess diagram for 
the RGB region of And XXI. The color represents the probability, Pcmd m each cell for a star in that position belonging to And XXI. 
This diagram is used as a probability map for all our spectroscopically observed stars. Red points represent stars DEIMOS stars for which 
Pcmd > 10 — 6 , while the blue stars represent objects for which Pcmd < 10~ 6 . 



Phalo 



exp 



\J 27r (^lhalo + V err,i) 



v r ,halo 




(4) 



halo 



The MW halo population is well approximated by 2 
Gaussians with v r< MWi = — 81.2 km s -1 , cr v ,MWi = 
36.5 km s -1 and v r ,MW2 = — 40.2 km s -1 and cr v ,MW2 = 
48.5 km s -1 , resulting in a probability density function 
for a given star with a velocity Vi and velocity uncer- 
tainty of v err i of: 



Pmw — 



R 



exp 




MWl 



V r ,MWl — v r,i 
V,MW1 V err,i / 

(1-R) 



(5) 



exp 



MWl 



1 / Vr,MW2 — v r,i 

2 



a v,MW2 + V err,i , 



where R is the fraction of stars in the first MW peak, and 
(1 — R) is the fraction of stars in the second peak. The 
value of R is determined empirically from our DEIMOS 
data set. 

A strong kinematic peak outside of these two popula- 
tions can then be searched for using a maximum likeli- 



hood technique, based on the approach of iMartin et al.l 
(2007). We search for the maximum in the likelihood 
function that incorporates the two contamination popu- 
lations plus an additional Gaussian structure with sys- 
temic velocity v r . su b s tr and a dispersion of (J v .substr, de- 
fined as: 



P. 



1 



substr 



exp 



^2it{[r}<7 VtSubstr \ 2 + v 2 errl ) 



ur, substr ^r,i 




(6) 



substr \ 



Here, to ensure we haven't biased our P su bstr strongly 
against stars that lie within the wings of the Gaussian 
distribution of dwarf spheroidal velocities, we have in- 
cluded a multiplicative free parameter, rj, to our derived 
value of <j v . To determine the ideal value of 77, we ran 
our algorithm over all our datasets, changing the value 
of 77 from 0.5 — 10.5 to see its effect on the final derived 
systemic velocities and velocity dispersion. We find that 
in all cases, the solutions converge at values of 77 ~ 2 — 4. 
For dSphs whose kinematics are well separated from con- 
taminants, the derived kinematics can remain stable up 
to much larger values of 77, however for those with sys- 
temic velocities within the velocity regime of the Milky 
Way, the solution quickly destabilizes as more contami- 
nants are included as probable members. We show this 
implicitly in Fig. 01 where we present the effect of modi- 
fying r) on 6 dSphs, And XII, XIX, XXI, XXII, XXIII 
and XXV. These objects were selected as they nicely 
probe our datasets with low numbers of probable mem- 
ber stars (~ 8), to those where we have 10s of probable 
members, as well as sampling dSphs from highly contam- 
inated to well isolated kinematic regimes. The value of 
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Fig. 3. — A velocity histogram for all observed stars in the field 
of And XXI. Our empirical Gaussian fits to the full Keck II data 
set for the M31 and MW halos are overlaid in blue and green 
respectively. A cold, kinematic peak at — 400 km s -1 is also seen, 
and this is the likely signature of the dSph. Our coarse, initial 
ML procedure identifies this peak, and the values of v syB and a v 
it measures are used to derive our kinematic probability, P ve i. 



rj is therefore independently determined for each dataset 
separately, and we report its final value in Tabled 



Psubstr 



1 



exp 




^^(^v.substr + V err,i) 



^r,substr ^r,i 



a 2 4- 7) 2 

v,substr 1 err,z 



(7) 



The likelihood function can then be simply written as: 



IN 

l0g[£(v r , <7 V )] = ^ l0 g( aP i,halo + /?Pi,MW + 7 P i, 



i.substr 



(8) 

where a, /3 and 7 represent the Bayesian priors, i.e., the 
expected fraction of stars to reside in each population. 
These are determined by starting with arbitrary fractions 
(for example, 0.2, 0.5 and 0.3 respectively) and are then 
adjusted to the posterior distribution until priors and 
posteriors match. This technique will therefore identify 
an additional kinematic peak, independent of the MW 
and M31 halo populations, if it exists. We stress that 
these are not the final systemic velocity and dispersion 
of the dSph, but merely indicate a region in velocity space 
in which an excess of stars above the two contaminant 
populations is seen. In Fig. [31 we show the result of this 
process for the And XXI dSph. Here, the substructure 
is clearly visible as a cold spike at ~ — 400 km s -1 . 

Now that we have a velocity profile for our three com- 
ponents (MW, M31 halo and the dSph), we can assign 
probabilities for each star within our sample belonging 
to each population using simple Bayesian techniques, i.e., 
the probability that a given star belongs to the substruc- 
ture, P„ei, is: 



aPfialo + PPmW + iPsubstr 

and the probability of being a contaminant is: 

_ aPhaio + /3Pmw 

aPhalo + PPMW + jPsubstr 



P, 



(9) 



(10) 



3.4. Measuring v r and o~ v 

Upon applying this to our data, we can identify the most 
probable members of each dSph, without having to apply 
any additional constraints or cuts. 

Having established the membership probability for 
each observed star (as detailed above) we now calculate 
the kinematic properties of each dSph; namely their sys- 
temic velocities (v r ) and velocity dispersi ons (a v ). We 
use th e maximum likelihood technique of iMartin et al.l 
(|2007D . modified to include our probability weights for 
each star. We sample a coarse grid in (u r , o~ v ) space and 
determining the parameter values that maximise the like- 
lihood function (ML), defined as: 

1 N _ 2 
]Dg[C(v r) a v )] = --J2\^og(a 2 tot )+P t ^-^ +P t \og(2n) 

2 f—f L otot 
i— 1 

(11) 

where N is the number of stars in the sample, v r .i is the 
radial velocity measured for the z th star, v err ^ is the cor- 

In this 



responding uncertainty and otot =- v < 

way, we aim to separate the intrinsic dispersion of a sys- 
tem from the dispersion introduced by the measurement 
uncertainties. 

3.5. Testing our probabilistic determination of 
membership and calculations of kinematic 
properties 

Having developed the above technique, it is important 
for us to rigorously test that it is robust enough to ac- 
curately determine the global kinematic properties for 
each of our datasets. In Appendix A, we examine in 
detail a number of potential issues that could cause our 
algorithm to return biased or incorrect results. These are 
the inclusion of a velocity dependent term in our prob- 
ability calculation, the effect of including low S:N data 

o 

(S:N< 5 A) in our analysis and the effect of small sample 
sizes (A* < 8) on our measurements of kinematic prop- 
erties. We briefly summarize our findings here, and refer 
the reader to Appendix A for a more detailed description. 

This work has introduced the concept of assigning a 
probability of membership for a given star to a dSph 
based on the prior knowledge of the velocity profiles of 
our expected contaminant populations, P ve i, a technique 
that has not previously been used in the study of M31 
dSphs. To test that this is not biasing our results, we 
can simply remove this term from Eqn. [TJ and follow the 
technique of T12, where they use only Pcmd and Pdist 
terms and then cut all stars with Pj < 0.1 and velocities 
that lie greater than 3<r from the mean of this sample 
from their final analysis. We find that both techniques 
produce very similar results, however our algorithm is 
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Fig. 4. — The effect of modifying r/ — a multiplicative weight applied to the dispersion of the identified substructure in our P ve i term - on 
the final derived systemic velocities (v r ) and dispersions (<j„) for six of our dSphs.The effect of increasing r] is typically more pronounced 
in objects with systemic velocities close to that of Milky Way contamination (such as And XIX and XXV) than for other objects. Dashed 
lines represent the optimal value of r) for each system. 



more robust in regimes where the systemic velocity of 
the dSph is close to that of the MW, and in dSphs where 
our number of probable member stars is low (A* < 10). 

To test the effect of low S:N data on our calculations 
of v r and a v , we use our datasets for And XXI, XXIII 
and XXV, all of which have > 25 associated members. 
For each dataset, we apply a series of cuts to the sample 

based on S:N (at levels of S:N> 2,3,4 and 5 A) and 
rerun our algorithm. In all cases we find that the derived 
probabilities do not significantly differ when the low S:N 
data are included, justifying our inclusion of all stars for 
which velocities are calculated by our pipeline. 

We also test the ability of our algorithm to measure v r 
and a v in the small regime. For some of our datasets, 



we are only able to identify a handful of stars as probable 
members. In theory, one can calculate velocity disper- 
sions accurately from only 3 stars if one is confident of 
ones measu r emen t uncertainties, as is demonstrated by 
lAaronsonTI (|1983f) measurement of the velocity disper- 
sion for Draco from only 3 stars, which remains consis- 
tent with moder n day measurements from significantly 
larger datasets ([Walker et all l2009b| ) . We can test if 
our results are similarly robust using our larger datasets 
(such as And XXI, XXIII and XXV) by randomly select- 
ing 4, 6, 8, 10, 15, 20 and 25 stars from these datasets 
and rerunning our algorithm to determine v r and <r v from 
these subsets. We repeat this exercise 1000 times for each 
sample size, and examine the mean and standard devia- 
tions for the computed quantities. We find that on aver- 
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age, for all sample sizes, our routine measures systemic 
velocities and velocity dispersions that are entirely con- 
sistent with those measured from the full sample, with 
a spread that is very comparable to typical errors pro- 
duced by our ML routine in these low A* regimes. As 
such, we conclude that our technique is able to place sen- 
sible limits on these values, even when dealing with as 
few as four member stars. 

Finally, as the individual positions, velocities and ve- 
locity uncertainties for all the stars analyzed in T12 are 
publicly available (with the data from the non-member 
stars having been kindly passed on to us by the SPLASH 
team) , we can check that our algorithm is able to repro- 
duce the values they measure for their M13 dSph sample. 
We find that, in all cases, we calculate systemic velocities 
and velocity dispersions that agree with their measured 
values to well within their la uncertainties. 

These tests demonstrate that our method is robust 
enough to accurately determine the global kinematic 
properties of M31 dSphs across a wide range of sample 
sizes and data quality. We therefore proceed to apply it 
to the datasets of all of the dSphs for which our group 
has acquired Keck II DEIMOS observations to date. 

4. THE KINEMATICS OF M31 DSPHS 

With this vast dataset of dSph kinematics now in hand, 
we can begin to statistically probe their structures more 
fully. In Fig. [5]we display a summary of the velocities and 
positions of 26 of the 28 dSphs for which kinematic data 
are avail able, where the value s are taken from this work, 
T12 and iKalirai et all (|2010f ) . In the following sections, 
we will discuss the individual stellar kinematics, masses 
and chemistries of the dSphs analyzed within this work. 

4.1. Andromeda XVII 

And XVII was discovered by llrwin et al.1 ((2008) , and 
it is located at a projected distance of ~ 40 kpc to the 
North West of Andromeda. A detailed study of deep 
imaging obtained with the Lar ge Binocular Camera o n 
the LBT was also performed by Brasseur et al.l (|2011bD . 
and throughout, we use the structural properties as de- 
termined from this work. It is a faint, compact galaxy 
(My = -8.61, r ha i f = 1.24' or 262^ pc). In the left 
panel of Fig. [5] we display the PAndAS color magnitude 
diagram for And XVII. Over-plotted we show the ob- 
served DEIMOS stars, color-coded by their probabilities 
of membership. The open symbols represent stars for 
which Pi < 10~ 6 . We employ this cut solely to make 
clarify which stars have the highest probability of be- 
longing to And XVII. In the right hand panel, we dis- 
play the basic kinematic information for And XVII. In 
the top panel of this subplot, we show a velocity his- 
togram for all stars observed within the LRIS mask, and 
stars with Pi > 10~ 6 are highlighted with a filled red his- 
togram. The centre panel shows the velocities as a func- 
tion of distance from the centre of And XVII (and the red 
dashed lines indicate 1,2,3 and 4 x n^if), and the lower 
panel shows the p hotometric metallic ities for all stars, as 
determined using iDotter et all (|2008| ) CFHT isochrones. 
Again, all points are color-coded by their probability of 
membership. Finally, the two lower panels show the re- 
sulting, one dimensional, probability weighted, marginal- 
ized maximum likelihood distributions for v r and a v for 
this data set. From the kinematics presented in Fig. [6l 



which represent the first spectroscopic observations of 
this object, we see the signature of the dwarf galaxy as a 
cold spike at v r ~ —250 km s^ 1 . From the lower panels of 
this figure and the accompanying CMD we see that there 
is a cluster of 7 stars sitting within this spike that are 
centrally concentrated and are consistent with the RGB 
of the dwarf itself, leading us to believe that our algo- 
rithm has cleanly detected the signature of the galaxy. 
Interestingly, we also see 3 stars that, kinematically, are 
indistinguishable from the stars that have been dubbed 
as probable members in our analysis. However, they all 
sit at large distances from the centre of And XVII, equiv- 
alent to greater than 6 times the half-light radius of the 
dwarf, and hence the routine has classified them as likely 
members of the M31 halo rather than And XVII mem- 
bers. But, given their tight correlation in velocity to the 
systemic velocity of And XVII the possibility exists that 
these are extra-tidal stars of And XVII. No sign of extra- 
tidal features were cited in either the discovery paper of 
And XVII or the LBT followup, but given its position in 
the north M31 halo where contamination from the MW 
becomes increasingly problematic, and its relatively low 
luminosity (My = —8.5), such features may be difficult 
to see within the imaging. However, at present they are 
considered unlikely members by our routine, and do not 
factor into our calculation of global kinematic proper- 
ties for this object. We find v r = — 251.6^2 kms -1 and 
a v = 2.9t^kms" 1 . 

4.2. Andromeda XVIII 

Andromeda XVI I I (A nd XVIII) was detected by 
iMcConnachie et al.l (|2008D in the PAndAS CFHT maps. 
Located at a projected distance of ~ 110 kpc to the 
North- West of M31, it is one of the most distant of its 
satellites, sitting ~ 600 kpc behind the galaxy, making 
spectroscopic observations of its individual RGB stars 
taxing, as they are all relatively faint (i > 22.2). Thus, 
from our 1 hour DEIMOS observation, we were only able 
to confirm 4 stars as probable members (see Fig. [7])- 
We determine the global systemic velocity to be v r = 
—346.8 ± 2.0kms _1 , and we are unable to resolve a ve- 
locity dispersion, finding a v = 0.0 +27 kms _1 where the 
upper bound is determined from the formal la confidence 
interval produced by our maximum likelihood analysis. 
This suggests that the 4 stars we are able to confirm as 
members do not adequately sample the underlying ve- 
locity profile. The systemic velocity we measure is dif- 
ferent to that presented in T12 of v r = —332.1 ± 2.7 at a 
level of 3.4cr. Our la limit of 2.7 km s -1 is also at odds 
with the dispersion determined by T12 {a v — 9.7 ± 2.3). 
They were able to measure velocities for significantly 
more probable member stars (22 vs. 4) owing to their 
longer integration of 3 hours. The faintness of our tar- 
gets and shorter exposure time could mean that strong 
sky absorption lines have systematically skewed our ve- 
locity measurements for these stars, and this could ex- 
plain the discrepancy of our measurements with respect 
to those of T12. To check against this, we again per- 
form weighted cross-correlations to each of the 3 Ca II 
lines individually. We find that the values we obtain, and 
their average, are fully consistent with that derived from 
the technique described in § 12.21 differing by less than 
3kms _1 from those values. The true systemic velocity 
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Fig , 5. — Positions of 26 of the 28 Andromeda dSphs, color coded by their heliocentric velocity, v r (taken from this work. iKalirai et alj 
120101 and T12). The solid line represents the PAndAS survey foot print. M31 is shown by the double ellipse, where the outer ellipse shows 
a segment of a 55 kpc radius ellipse flattened to c/a = 0.6, and the major-axis and minor-axis are indicated with straight lines out to this 
ellipse. The inner ellipse corresponds to a disc of radius 2° (27 kpc), with the same inclination as the main M31 disc. M33 is shown as a 
singular ellipse. 



of And XVIII therefore remains unclear. However, given 
their larger sample size, the T12 systemic properties are 
more statistically robust than those we present here. 

4.3. Andromeda XIX 

Andromeda XIX (And XIX) was first reported in 
IMcConnachie et"aLl (|2008D . and is a relatively bright, 
very extended (My = —9.3, r^aif = 1.5 kpc) dSph, 
located at a projected distance of ~ 180 kpc to the 
south west of M31. Its unusual morphology, very low 
surface brightness £„ = 30.2mag/arcsec 2 , and evidence 
in the photometry for a p ossible link to the major 
axis substructure re ported in llbata et al.1 f|2007T ) caused 
IMcConnachie et all (|2008l) to question whether And XIX 
was truly a dynamically relaxed system, or whether it 
had experienced a significant tidal interaction. Here, we 
present the first spectroscopic observations of the And 
XIX satellite in Fig. [5] from two DEIMOS masks placed 
at different position angles. These data allow us to com- 
ment on its dark matter content, and on the likelihood 
of a tidal origin for its unique structure. We identify 
27 stars where Pi > 10~ 6 within the system. These 
measurements were made increasingly challenging as the 



systemic velocity we measure is v r = — 111. 6^' ^kms -1 , 
placing it within the regime of Galactic contamination. 
However, we are confident that our algorithm is robust 
to this unfortunate location of And XIX in velocity space 
(see discussion in § 13.51 and Appendix A) . As a further 
check that none of the stars we define as probable mem- 
bers are actually foreground contaminants, we measure 
the equivalent widths of the Na I doublet lines (located 

o 

at ~ 8100 A). These gravity-sensitive absorption lines 
are typically significantly stronger in foreground dwarf 
stars than M31 RGB stars, although there is some over- 
lap between the two populations. For the stars tagged 
as probable members by our algorithm, we find no ev- 
idence of strong absorption in the region of the Na I 
doublet, indicating that we are not selecting foreground 
stars as members. We measure a relatively cold velocity 
dispersion for this object of a v = 4.71}' 4 kms -1 , which is 
surprising given the radial extent of this galaxy. This re- 
sult will be discussed further in § 15.11 and in a follow-up 
paper (Collins et al. in prep). 

4.4. Andromeda XX 
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Fig. 6. — CMD of And XVII, from PAndAS photometry (not extinction corrected). All stars observed with DEIMOS are color coded by 
their probability of membership. Top right: Kinematic information for all stars observed with DEIMOS. The top subplot shows a velocity 
histogram for our sample, with all stars with Pj > 10~ 6 highlighted as a filled red histogram. The central subplot shows the distance of 
each star from th e centre of And XV II as a function of velocity, and the lower subplot shows the photometric metallicity for each star, 
interpolated from Dottcr ct al. (2008) i sochro nes in the CFHT-MegaCam g— and i— bands (after correcting the observed colors of the stars 
for extinction using the ISchlegel et al.1 $1993 ) dust maps), as a function of velocity. Stars that lie far from the RGB of the dSph are not 
well matched by the isochrones, and as a result there is no estimate of their [Fe/H]. Lower panels: Resulting probability distributions 
from ML analysis of the kinematics of And XVII. The left hand panel shows the systemic velocity (v r ) likelihood distribution and the right 
shows the intrinsic velocity dispersion (cr v ) likelihood distribution. The dashed lines represent canonical 1,2 and 3<r confidence intervals, 
derived assuming Gaussian uncertainties. 
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And XX was the third of three dSphs discovered by 
iMcConnachie etldl (|2008f) . and is notable for being one 
of the faintest dSph companions detected surrounding 
Andromeda thus far. With My — —6.3 and n ia if = 
114^2 P c : ^ is a challenging object to study spectro- 
scopically as there are very few stars available to tar- 
get on its RGB, as shown in the top left subplot in 
Fig. [9l As a result, our algorithm is only able to find 
4 stars for which P, > 10 -6 . These are found to clus- 
ter around v r = — 456.2l 3 ; 6 kms -1 , with a dispersion of 



a v = 7. 1+25 kms -1 . Despite the low number of stars, we 
are confident in this detection, as the systemic velocity 
places it in the outer wings of the velocity profile of the 
M31 halo. And XX is also located at a large projected 
distance from M31 of ~ 130 kpc, where we expect the 
density of the M31 halo to be very low. As such, seeing 4 
halo stars so tightly correlated in velocity in the wings of 
the halo velocity profile within such a small area of the 
sky (all stars are within 1 arcmin of the centre of And 
XX) is highly unlikely. We caution the reader that, while 
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we are confident that our algorithm is able to measure 
velocity dispersions for sample sizes as small as 4 stars, 
as we are not probing the full velocity profile of this ob- 
ject this measurement ideally needs to be confirmed with 
larger numbers of member stars. 

4.5. Andromeda XXI 

Andromeda XXI (And X XI) was identified wit hin the 
PAndAS imaging maps by IMartin et all (|2009f ). It is 
a relatively bright dSph (My = —9.9), located at a 
projected distance of <~ 150 kpc from M31, and it 
has a half-light radius of ri ia if = 842 ± 77 pc. We 
present our spectroscopic observations for this object in 
Fig. [TUl and in the top right subplot, we can clearly 
see the signature of And XXI as a cold spike in ve- 
locity with 32 probable member stars, located at v r = 
—362.5 ± 0.9 km s -1 , with a curiously low velocity dis- 
persion of only <r v = 4.5^'q kms -1 . These results are 
completely consistent with those of T12, where they mea- 
sured v r — — 361.4±5.8 km s —1 and o~ v = 7.2±5.5kms~ 1 . 
As their sample contained only 6 likely members com- 
pared with the 29 we identify here, ours constitute a more 
statistically robust measurement of the global kinematics 
for this object than those presented in T12. 

4.6. Andromeda XXII 

Andromeda XXII (And X XII) was ident i fied w ithin the 
PAndAS imaging maps bv IMartin et al.l fj2009f ) . and is 
a relatively faint dSph, with My = —6.5. Its physical 
position in the halo, located at a distance of 224 kpc 
in projection from M31, but only 42 kpc in projection 
from M33, led the authors to postulate that it could be 
the first known dSph satellite of M33. Subsequent work 
analysing the kinematics of And XXII by T12 measured 
a systemic velocity for And XXII of —126.8 ± 3.1 kms -1 
from 7 stars, more com patible with t he systemic velocity 
of M33 f-178k ms- 1 . lMateolH998l) tha n that of M31. 
Another study bv lChapman et al.l ([2012D using the same 
data and the same method we present here concluded the 
same, measuring a systemic velocity for the satellite of 
— 129.8 ± 2.0 from 12 probable member st ars, consistent 
with the T12 value. IChapman et al.l (|2012l ) also compare 
the position and kinematics of And XXII with a suite of 
N— body simulations of the M31-M33 system, concluding 
that And XXII was a probable M33 satellite. 

I n Fig. ITTl we present the same data as analyzed 
by IChapman et al.l <|2012f ) for completeness. The ve- 
locity dispersion of And XXII is just resolved at a v = 
2.8l 1 '4kms~ 1 , completely consistent with the value of 
a v = 3.5^2'1^ ms l from T12. As our values are calcu- 
lated from a 50% greater sample size, we posit that they 
are the more statistically robust. 

4.7. Andromeda XXIII 

Andromeda XXIII (An d XXIII) was the first of five 
M31 dSphs identified bv iRichardson et al.l (j20Fu). Lo- 
cated at a projected distance of ~ 130 kpc to the east 
of Andromeda, it is relatively bright, with My — —10.2, 
and extended, with rhaif = lOOlt^ P c - Our routine 
clearly detects a strong cold kinematic peak for And 
XXIII located around —230 kms -1 and calculates a sys- 
temic velocity of v r — —237.7 ± 1.2 kms -1 , and a veloc- 
ity dispersion of a v — 7.1 ± 1.0 kms -1 from 40 probable 



member stars, as show in Fig. 1121 This small, positive ve- 
locity relative to M31, combined with its large projected 
distance from the host suggests that And XXIII is not 
far past the apocentre of its orbit, heading back towards 
M31. 

4.8. Andromeda XXIV 

And XXIV was also first reported in IRichardson et "al~l 
(|2011t ). Relatively faint and compact (My = —7.6, 
Hiaif = 548^37 pc), spatially it is located ~ 200 kpc from 
M31, along its northern major axis. And XXIV was ob- 
served on two separate occasions as detailed in Table [5J 
For the first mask, there was an error in target selection, 
and as a result, only one star that lay on the RGB of And 
XXIV was observed. The second mask was observed in 
May 2011, however owing to target visibility, only a short 
integration of 45 minutes was obtained, which resulted 
in higher velocity uncertainties than typically expected 
(~ 8kms _1 vs. ~ 5kms -1 ). For this reason, we have 
only included stars from this mask with i < 22.0, as the 
spectra for fainter stars were too noisy to determine reli- 
able velocities from. The systemic velocity of And XXIV 
also unsatisfactorily coincides with that of the MW halo 
contamination, as can be seen in Fig. [13l As for And 
XIX, we check the strength of the Na I doublet of all 
the stars classified as potential members for And XXIV, 
and find no significant absorption, making them unlikely 
foreground contaminants. But, given the lower quality 
of this dataset, this check is far from perfect, and it is 
possible that we have included contaminants from the 
MW within our sample. Owing to the larger velocity un- 
certainties of the And XXIV dataset, and the overlap of 
And XXIV with the MW, the determination of probabil- 
ity of membership for stars within this dataset is based 
largely on their position in the color magnitude diagram 
(e.g., location on the RGB) and their distance from the 
centre of And XXIV. 

When we run our machinery over the data acquired 
from both masks, we identify only 3 probable members 
and determine a systemic velocity of v r = —128.2 ± 
5.2 km s -1 and we resolve a velocity dispersion of <r v — 
0.0 +7,3 kms" 1 . Given the lower quality of this dataset in 
comparison to the remainder of those we present in this 
work, and the overlap of And XXIV in velocity space 
with contamination from the MW, a robust kinematic 
detection and characterisation of this galaxy is made in- 
credibly challenging. As such, we present these results 
as a tentative identification of And XXIV, and do not 
include its measured properties in the remainder of our 
analysis. Further kinematic follow up of And XXIV is re- 
quired to understand this system. We present the veloc- 
ities of all bright stars for which velocity measurements 
were possible in Table [3] so that they may be helpful for 
any future kinematic analysis of this system. 

4.9. Andromeda XXV 

And XXV was identified in IRichardson et al.1 (|2011f ) 
as a relatively bright (My = —9.7), extended (r^aif = 
642^74 pc) dwarf spheroidal, located at a projected dis- 
tance of ~ 90 kpc to the north west of M31. As with And 
XXIII, we present here a kinematic analysis of And XXV. 
The results are displayed in Fig. [14] We see that the sys- 
temic velocity of And XXV (v r = -107.8 ± l.Okms" 1 ), 
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Fig. 9.— As Fig.© but for And XX. 



places it in the regime of the Galactic foreground. How- 
ever, given the strong over-density of stars with this ve- 
locity relative to the expected contribution of MW stars, 
we are confident that our routine has detected 25 likely 
members for this object. We check the strength of the Na 
I doublet in these likely members, and find no significant 
absorption, making them unlikely foreground contami- 
nants. As for And XIX and XXI, we find that And XXV 
has a curiously low velocity dispersion for its size, with 
<j v = 3.0i^ '^kms -1 . We discuss the significance of this 



further in § 15.11 and Collins et al (2013, in prep). 
4.10. Andromeda XXVI 
And XXVI is a relatively faint (My = -7.1) dSphwith 
Hialf = 219^52 pc, also first reported in lRichardson et aTl 
(|2011| ). Its low luminosity makes observing large num- 
bers of member stars difficult, owing to the paucity 
of viable targets on the RGB that can be observed 
with DEIMOS. As a result, our routine has identi- 
fied only 6 stars as potential members, highlighted in 
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Fig. 10.— As Fig. [1 but for And XXI 



Fig. [15] The dwarf has a systemic velocity of v r = 
— 261.6i|" gkms -1 , and a fairly typical velocity disper- 
sion of a v = 8.6i 2 '2 kins" 1 . As with And XX, while we 
believe our routine can robustly measure the velocity dis- 
persions of systems with only 6 confirmed members, to 
be truly confident of this value, follow up of And XXVI 
to i ncrease the number of likely members is required. 

In iConn et al.l (|2012f ) , from an analysis of the photom- 
etry of And XXVI, they determined a distance modulus 



to the object of (m - M) = 24.39t° j,l from a Markov- 
Chain-Monte-Carlo analysis of the PAndAS photometry 
of And XXVI. This value corresponds to an i— band mag- 
nitude for the TRGB of And XXVI of m ifi = 21.llg;||. 
Our CMDs for the dwarfs are not ext inction corrected, 
but u sing the extinction valu es from iRichardson et al.l 
(|20T1 of E(B - V) = 0.110 (|Schlegel et al II1998T) . this 
would correspond to an i— band magnitude of rrii : TRGB = 
21.3t£ff. Three targets were observed with magnitudes 
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Fig. 11.— As Fig. [6] but for And XXII 



and colors that should be consistent with their belonging 
to And XXVI. However, we find that all these objects 
have velocities that are consistent with being Galactic 
foreground contaminants. Given the position of And 
XXVI in the northern M31 halo, where contamination 
from the MW increases, this is not unexpected. The 
brightest star we observe that is likely associated with 
And XXVI has mi = 21.9 (m ifl = 21.7). Assuming that 
this brightest confirmed member star of And XXVI sits 
at the RGB tip, the distance estimate becomes 1.1 Mpc. 



While this value is higher than that of lConn et al.l (|2011f) , 
it is still within their upper distance estimate. Addition- 
ally, as not every star brighter than the member with 
nii.o = 21.7, with colors consistent with the And XXVI 
RGB was observed, this value merely represents an upper 
limit, on the distance to And XXVI and highlights the 
difficulty of calculating distances to these faint galaxies 
where RGB stars are sparse. 

4.11. Andromeda XXVII 
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Fig. 12.— As Fig. [1 but for And XXIII. 

Andromeda XXVII (And XXVII) is a somewhat un- 
usual object as it is currently undergoing tidal disrup- 
tion, spreading its constituent stars into a large stellar 
stream, named the northwestern arc , disco vered in the 
PAndAS survey by Richard son et al.1 (|2011l ). As such, it 
is unlikely to be in virial equilibrium, if it remains bound 
at all. 

When determining the kinematics of And XXVII, we 
find the results somewhat unsatisfactory. Our rou- 
tine determines v r = — 539. 6^45 kms -1 and a v — 



5 10 15 20 25 30 
a v (kms 1 ) 



14.8^3'^ kms -1 from 11 stars. However, from an inspec- 
tion of Fig. [TBI we see that there is significant substruc- 
ture around v hei ~ —500 km s _1 , much of which is consid- 
ered to be unassociated with And XXVII in this analysis 
as it does not fall within a cold, well-defined Gaussian ve- 
locity peak. Given the disrupting nature of And XXVII, 
it is likely that a different analysis is required for this ob- 
ject, and we shall discuss this further in a future analysis, 
where the kinematics of the northwestern arc itself are 
also addressed. From this first pass however, it would ap- 
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TABLE 3 

Details of probable dSph members (This table is available in its entirety in machine-readable form in 
the online journal. a portion is shown here for guidance regarding its form and content) 



Field 


Star ID 


"0,J2000 


<5o,J2000 


a 


i 


v hel (kms 1 ) 


» err (kms 1 ) 


S:N ( A" 1 ) 


Pi 


And V 


9 


1:10:2.38 


47:37:48.5 


23.640 


22.402 


-371.270 


5.080 


1.700 


0.014 


And V 


12 


1:10:5.540 


47:36:41.6 


23.100 


21.530 


-406.500 


3.600 


2.000 


0.111 



pear that And XXVII may no longer be a gravitationally 
bound system. 

4.12. Andromeda XXVIII 

And XXVIII was recently discovered in the 8th data 
release of the SDSS survey ([Slater et all 120111 ). It has 
M v = —8.5 and 7*h a if = 210^5° P c - K is also potentially 
one of Andromeda's most distant satellites, with a host- 
satellite projected separation of 365^ 7 kpc. The And 
XXVIII satellite is not covered by the PAndAS footprint, 
so we must instead use the original SDSS photometry for 
our analysis. A CMD with the SDSS i— band and r — i 
colors for And XXVIII is shown in Fig. [T71 where all 
targets brighter that i ~ 23.5 within 3rh a if are shown. 
The photometry here do not show an RGB that is as 
convincing as those from the PAndAS surve y, so to guide 
the ey e, we also overplot an isochrone from [Potter et afl 
(20081) with a metallicity of [Fe/H] = -2.0, corrected for 
the di stance of And XXVIII as reported in iSlater et al.1 
([20T1 . 

In a recent paper, iTollerud et al.l ([2013) discussed the 
kinematics of this object as derived from 18 members 
stars. They find v r = — 328.0 ± 2.3 kms -1 and a v = 
8. 1 ± 1.8 km s ~ 1 from their full sample . They then remove 
two stars that they categorize as outliers based on their 
distance from the centre of And XXVIII, which alters 
their measurements to v r = — 331.1±1.8kms and o~ v — 
4.9 ± 1.6kms -1 . Analyzing our own PEIMOS dataset 
for this object, we find v r = —326.1 ± 2.7 km s^ 1 and 
a v = 6.6^2 2 s^ 1 based on 17 probable members. This 
is fully consistent with the results from the full sample 
in ITollerud et al.l (|2013|) . However, the systemic velocity 
we measure is offset at a level of ~ la from their final 
value (calculated after excluding 2 outliers). This offset 
is small, and is probably attributable to our differing 
methodologies for classifying stars as members. As we 
believe our method is more robust (as discussed in § 13.51 
and Appendix A), we will use our derived parameters for 
this object in the remainder of our analysis. 

4.13. And XXX/Cassiopeia II 

And XXX - also known as Cass II owing to its spatial 
location, overlapping the Cassiopeia constellation - is a 
recently discovered dSph from the PAndAS survey (Irwin 
et al. in prep). It has M v = —8.0 and rhaif = 267lgg pc. 
Located to the north west of Andromeda, it sits within 
60 kpc of the two close dwarf elliptical M31 companions, 
NGC 147 and NGC 185. With these 3 objects found so 
close together in physical space, it is tempting to suppose 
them a bound system within their own right, but this can 
only be borne out by comparing their kinematics. 

Conspiring to confound us, we find that Cass II has 
kinematics that place it well within the regime of Galac- 



tic foreground, as can be seen in Fig. [T5] However, 
our analysis is able to detect the dSph as a cold spike 
consisting of 8 likely members. As for And XIX, we 
check the strength of the Na I doublet in these likely 
members, and find no significant absorption, making 
them unlikely foreground contaminants. We measure 
v r = — 139.8lgg kms -1 , and a fairly typical velocity dis- 
persion of a v — 11. kms -1 . 

The systemic velocity of Cass II (v r = 
-lSg.Stg^kms- 1 ) puts it within ~ SOkms" 1 of 
those of NGC 147 and NGC 185 (v r = -193 ± 3 kms" 1 
and v r = -210 ± 7kms _1 , IMateol Il998f ). lending 
further credence to the notion that these 3 systems are 
associated with one another. This will be discussed in 
more detail in Irwin et al. (2013, in prep). 

4.14. A note on previous work 

Finally, we also use our new algorithm to reanalyze all 
our previously pu blished M31 dSph d atasets. These in- 
clude And V. VI (ICollins et all [2011 . X I, XII and XIII 
(|Chapman et al.l2007t ICollins et al.ll2010D . Petailsofthe 
results of this reanalysis can be found in Appendix B. In 
summary, we find that our algorithm measures systemic 
velocities and velocity dispersions that are fully consis- 
tent with our previous work. We present these results in 
Table [U And V, XI, XII and XIII are also analyzed by 
T12, so we compare our findings with theirs. For And 
XI and XII, our results are based on two and four times 
the number of stars respectively, and as such, supercede 
those presented in T12. In the case of And V and XIII, 
the T12 measurements are based sample sizes with four 
times the number of stars as our datasets, making their 
findings more robust. 

In previous stud i es by our group (iChapman et al.l 
120051: ILetarte et alj [2009t ICollins et all 12010ft . we also 
published kinematic analyses for three additional M31 
dSphs; And IX, And XV and And XVI. In T12, it was 
noted that the values presented in these works for sys- 
temic velocities and velocity dispersions were not con- 
sistent with those measured in their analyses. We re- 
visited these datasets in light of this discrepancy, to see 
if our new technique could resolve this issue. We found 
that these discrepancies remained. For And IX, we mea- 
sure a systemic velocity of v r — —204.8 ± 2.1 kms -1 cf. 
v r = — 209. 4±2. 5 km s _1 in T12 and a velocity dispersion 
of a v = 2.0t^kms" 1 cf. a v = 10.9 ± 2.0kms- 1 . Not 
only are their measurements determined from 4 times the 
member stars that we possess, we also experienced prob- 
lems with our radial velocity measurements for the stars 
observed with this mas k, due to the use of the minislitlct 
approach pioneered by Ib ata et all (|2005f h This setup 
resulted in poor sky subtraction for many of the science 
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Fig. 13.— As Fig. [1 but for And XXIV. 



spectra, lowering the quality of our radial velocity mea- 
surements. As such , the T12 results supe r cede those of 
our p revious work (|Chapman et al.l 120051 : iCollins et al.l 
I201TO . 

For And XV and XVI, we measure a systemic velocities 
ofv r = -354.6i4.9kms" 1 and v r = -374.1±6.8 kms" 1 
cf. v r = -323.0 ± 1.4 kms" 1 and v r = -367.3 ± 
2.8 kms" 1 from T12. We also note offsets in our ve- 
locity dispersions for And XV and XVI, where we mea- 
sure a v = 9.61^6 kms -1 and a v — 17. 3^ 4 kms" 1 cf. 



<j v = 3.8 ±2.9 kms" 1 . In this instance, the data for both 
And XV and XVI were taken in poor conditions, with 
variable seeing that averaged at 1.8" and patchy cirrus. 
These conditions significantly deteriorated the quality of 
our spectra, and made the measurement of reliable ra- 
dial velocities extremely difficult. Again, this leads us to 
conclude that the measurements made in T12 superc ede 
those presented by our group in lLetarte et al.l (|2009D . 
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Fig. 14.— As Fig. [1 but for And XXV. 

5. THE MASSES AND DARK MATTER CONTENT OF M31 

DSPHS 

5.1. Measuring the masses and mass-to-light ratios of 
our sample 

As dSph galaxies are predominantly dispersion sup- 
ported systems, we can use their internal velocity dis- 
persions to measure masses for these systems, allowing us 
to infer how dark matter dominated they are . There are 
several methods in the literature for this (e.g. Jlllingworthl 



5 10 15 20 25 
a v (kms 1 ) 



| 1976t l Richsto ne"fc Tremainelll986l) . but recent work by 
IWalker et alT (f2009a) has shown that the mass contained 
within the half-light radius (Mh a if)of these objects can 
be reliably estimated using the following formula: 

Mhalf = MHialfC^half ( 12 ) 

where /i = 580 M pc _1 km _2 s 2 , rtaif is the spherical half- 
light radius in pc and avhaif is the luminosity- averaged 
velocity dispersion. This mass estimator is independent 
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Fig. 15.— As Fig. El but for And XXVI. 



of the (unknown) velocity anisotropy of the tracer pop- 
ulation, however, it is sensitive to the embeddedness of 
the stellar component within the DM halo. Particularly, 
the mass tends to be s lightly over-estimated the m ore 
embedded the stars are (Walker & Pcharrubia 2011), es- 
pecially if the dark matter halo follows a cored density 
profile. 

As numerous authors have shown that the veloc- 
ity dispersio n profiles of dSph s are c onstant with 
radius (e.g., IWalker et all 120071 l2009bf ). we assume 



our measured values of <j v are representative of the 
luminosity- averaged velocity dispersion (cr^haif) used by 
IWalker et al.l (|2009bD . However, if it transpired that the 
velocity dispersion profiles of the Andromedean dSphs 
were not flat, but declined or increased with radius, this 
would no longer true. We see no evidence for this be- 
haviour in our dataset, although low-number statistics 
means we are unable to completely rule out this pos- 
sibility. We calculate this for all our observed dSphs 
(including those we reanalyzed from previous works, see 
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Fig. 16.— As Fig. [1 but for And XXVII. 



Appendix B) using results from the Keck LRIS and 
DEIMOS dataset, and report their masses within r na if 
(M ha i f ) in Table H 

From these masses, it is trivial to estimate the dynami- 
cal central mass-to- light ratios for the objects, [M/L]h a if • 
We list these values for each dSph in Table IH where the 
associated uncertainties also take into account those from 
the measured luminosities and distanc es to these dSphs 
(|McConnachiell2012t iConn et all 12012ft . as well as those 
on the masses measured in this work. 



5.2. Comparing the mass-to-light ratios of MSI and 
Milky Way dSphs 

By combining our measurements of the kinematic 
of M31 dSphs in thi s work with those from T12 and 
iTollerud et all p013ft . we find ourselves with a set of 
kinematic properties as measured for 27 of the 28 An- 
dromeda dSphs (owing to the difficulties experienced 
with the And XXIV dataset, we do not include this 
object in our subsequent analysis). This near-complete 
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Fig. 17.— As Fig. [1 but for And XXVIII. 



sample allows us to fully compare the masses and mass- 
to- light ratios for the M31 satellite system with those 
measured in the Milky Way satellites. Before beginning 
this analysis, we compile Table[5]which presents the kine- 
matics for the full M31 satellit e system, which com bines 
the results from this work, T12. lKalirai et all (|2010f ). and 
iTollerud et all (|2013D . In cases where two measurements 
for a dSph exist, we use those that were calculated from 
larger numbers of likely members, as these are the more 
robust. We begin by comparing the mass-to-light ra- 



tios (which indicate the relative dark matter dominance 
of these objects) of the two populations as a function 
of luminosity. In Fig [19] we show these values for all 
MW ( red triangles, with values taken from I Walker et akl 
I20r)9bn . and M31 (blue circles) dSphs as a function of 
their luminosity. We can see that all these objects are 
clearly dark matter dominated, excluding And XII and 
And XVII where we are unable to resolve the mass with 
current datasets. We also see that they follow the trend 
of increasing [M/L] na if with decreasing luminosity, as is 
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Fig. 18.— As Fig. [6] but for Cass II 



seen in their MW counterparts. 

The one tentative exception to this is the And XXV 
dSph. From our dataset, we measure a value of 
[M/L]hait = 10.3lg7 f° r this object, making it consis- 
tent with a stellar population with no dark matter within 
its la uncertainties. This result is surprising and would 
be of enormous importance if confirmed with a larger 
dataset than our catalog of 26 likely members as it would 
be the first dSph to be observed with a negligible dark 
matter component. And XXV is also one of the mem- 



bers of the r ecently discovered thin plane of satellites in 
Andromeda (|Ibata et al.ll2013f) . and so the presence or 
absence of dark matter in And XXV might tell us more 
about the origins of this plane which are currently poorly 
understood. 

5.3. Comparing the masses of M31 and Milky Way 

dSphs 

Finally, we discuss how the masses for the full sam- 
ple of Andromeda dSphs for which kinematic data are 
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TABLE 4 

Kinematic properties of Andromeda dSph galaxies as derived within this work, 
from Keck I LRIS and Keck II DEIMOS data. 
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Fig. 19. — Dynamical mass-to-light ratio within the half-light 
radius, [M/LJ^aif, as a function of half-light radius for all M31 (blue 
circles), MW (red triangles) and isolated dSphs (cyan squares). 

available compare with those of the MW dSphs. For the 
M31 dSph population, we again use our compilation of 
kinematic properties assembled in Table [S] We plot the 
velocity dispersions, mass within the half-light radius, 
and central den sities for all M31 (blue circles) and MW 
(red triangles IWalker et all 120091* lAden et al.l 120091: 
iKoposov et alj|201lHSimon et al. ||2011D dSphs" as a func- 
tion of radius. We then overplot the best- fit NFW and 
cored mass profiles for the MW, taken from IWalker et aLI 



(|2009bD . In general, we see that the M31 and MW 
objects are similarly consistent with these profiles, an 
agreement that was also noted by T12. However, there 
are 3 objects which are clear outliers to these relations. 
These are And XIX, XXI and XXV, with velocity dis- 
persions of a v = 4.7+I 4 kms -1 , a v = 4.511'q kms^ 1 and 
a v = 3. Oil 1 kms -1 , as derived in this work. Given their 
half-light radii (r ha i f =1481±g 8 pc, r half =842±77 pc 
and rhaif =6421^4 pc), one would expect them to have 
dispersions of closer to 9 km s _1 in order to be consistent 
with the MW mass profile. As they stand, these three 
objects are outliers at a statistical significance of 2.5a, 
3.0a and 3.4cr (calculated directly from their likelihood 
distributions as presented in Figs. HI [TU] and [U|). Sim- 
ilarly, in T12 they noted that And XXII and And XIV 
were outliers in the same respect as And XIX, XXI and 
XXV, albeit at a lower significance. These difference can 
also be observed in terms of the enclosed masses and 
dens ities within H i aif- 

In iCollins etldl (|201l1) we argued that the low veloc- 
ity dispersion seen in some Andromeda dwarfs were a 
result of tidal forces exerted on their halos by the host 
over the course of their evolution, and that this effect 
was predominantly seen in dSphs where their half-light 
radii were more extended for a given luminosity than ex- 
pected, such is the case our three outliers, And XIX, XXI 
and XXV. This result therefore adds weight to the trend 
presented in that work. A number of recent works trying 
to account for the lower than predicted central masses of 
dSph galaxies within t he Local Group a ls o sup port this 
notion. For example, iPenarrubia et all (|2010T > demon- 
strated that the presence of a massive stellar disk in the 
host galaxy (such as those of the MW and M31) can sig- 
nificantly reduce the total masses of its associated satel- 
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Fig. 20. — Top: Half-light radius (r^aif) vs. velocity dispersion for the MW (red triangles) and M31 (blue circles) dSphs. The best 
fit NFW and cored mass profiles to the whole Local Group dSph population are over-plotted as dashed and dot-dashed lines respectively. 
Middle: As above, but in the Mhaif — Hialf plane. Lower: As above but in the phalf — ''half plane. 



lite s. In addition, recent, pap ers bv lZolotov eTaLl (pOll ) 
and lBrooks fc Zolotovl (|2012f ). where the effect of baryons 
within dark matter only simulations was measured also 
find that tidal forces exerted by host galaxies where a 
massive disk is present will serve to reduce the masses of 
its satellite population at a far greater rate than hosts 
without baryons. And XIX, XXI and XXV may thus 
represent a population of dSph satellites whose orbital 
histories about M31 have resulted in substantial fractions 
of their central mass being removed by tides. It should 
be noted, however, that tides not only reduce the central 
masses and densities of dSph halos, th ey also reduce the 
spatia l size o f the luminous component ([Penarrubia et al.l 
2008b, 2010f), albeit at a slower rate. The tidal scenario 
is therefore slightly difficult to reconcile with these outly- 
ing M31 dSphs having the largest sizes, unless they were 



both more massive and spatially larger in the past. 

Other recent theoretical works have also shown that 
the removal of baryons from the very centres of dark 
matter halos by baryonic feedback (from star forma- 
tion and supernovae, for example) can also help to 
lower the central masses and densiti es of satellite galax- 
ies fe.g.JPontzen fc Goyernatol[20il iZolotov et al.ll2012t 
iBrooks fc Zolotovl l2012f T" For this method to work ef- 
fectively, however, very large 'blow outs' of gas are re- 
quired, of the order ~ 10 8 — 1O 9 M0, equivalent to 
~ 40000 SNe. This would require a minimum initial 
satellite luminosity of My < -12 (|Zolotov et al.l 120121: 
iGarrison-Kimmel et al.ll2013f l , significantly brighter than 
the current luminosities of our outliers (My ~ —10). 
Therefore, if feedback has indeed played a role in the 
shaping of the dark matter halos of And XIX, XXI and 
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XXV, one assumes it would have to have operated in 
tandem with tidal stripping. We will discuss the im- 
plications and interpretation of these result further in a 
companion paper (Collins et al. in prep). 

6. METALLICITIES 

Our observational setup was such that we cover the 
calcium triplet region (Ca II) of all our observed stars. 
This strong, absorption feature is useful not only for 
calculating velocities for each star, but also metallici- 
ties. For RGB stars, such as we have observed, there 
is a well known relation between the equivalent widths 
(EWs) of the Ca II lines, and the iron abundance, 
[Fe/H], of the object. The calibration between these 
two values has been studied and tested by numerous au- 
thors, using both globular clusters and dSphs, and is 
valid down to metallicitie s as low as [Fe/H] ~ — 4 (see 
e.g- lBattadia et al.ll2008HStarkenburg et alM OlO). Fol- 
lowing the iStarkenburg et al.l ()201Cll ) method, which ex- 
tends the sensitivity of this method down to as low as 
[Fe/H] ~ —4, we fit Gaussian functions to the three Ca 
II peaks to estimate their equivalent widths (EWs), and 
calculate [Fe/H] using equation [T3l 

[Fe/H] = -2.87 + 0.195(Vr GB - V HB ) + 0.48ECa 

-0.913ECa~ 15 + 0.0155SCa(V RGB - V HB ) 

(13) 

where ECa=0.5EW 8 498+EW 85 42+0.6EW 8662 , V RGB is 
the magnitude (or, if using a composite spectrum, the 
average, S:N weighted magnitude) of the RGB star, and 
Vhb is the mean V- magnitude of the horizontal branch 
(HB). Using Vhb — Vrgb removes any strong depen- 
dence on distance or reddening in the calculated value of 
[Fc/H], and gives the Ca II line strength at the level of 
the HB. For M31, w e set this value to be Vhb =25.17 
(jHolland et al.lll996lPI . As the dSphs do not all sit at 
the same distance as M31, assuming this introduces a 
small error into our calculations, but it is at a far lower 
significance than the dominant uncertainty introduced 
by the noise within the spectra themselves. For individ- 
ual stars, these measurements carry large uncertainties 
(> 0.4 dex), but these are significantly reduced when 
stacking the spectra into a composite in order to mea- 
sure an average metallicity for a given population. 

Uncertainties on the individual measurements of 
[Fe/H] from our stellar spectra are typically large (> 
0.5 dex), so for a more robust determination of the av- 
erage metallicities we co-add the spectra for each dSph 

(weighting by the S:N of each individual stellar Spee- 
ch 

trum, which is required to be a minimum of 2.5 A -1 ) 
and measure the resulting EWs. In a few cases, not all 
3 Ca II lines are well resolved. For And V, IX, XVII, 
XVIII XXVI and XXVIII, the third Ca II line is signif- 
icantly affected by skylines, whilst for And XXIV, the 
first Ca II line is distorted. In the case of And XIII, 
only the second Ca II line appears well resolved. In 
these cases, we neglect the affected lines in our esti- 
mate of [Fe/H], and derive reduced equivalent widths 

1 5 This assumed v alue is sensitive to age and metallicity effects, 
sec Chen ct al. 2009 for a discussion, however owing to the large 
distance of M31, small differences in this value within the M31 
system will have a negligible effect on metallicity calculations 



from the unaffected lines. Where the third line is af- 
fected, this gives SCa=1.5EW 8 498+EW 8 542- Where the 
second line is affected, we find ECa=EW 8 542+EW 8 662- 
Finally, where only the second line seems reliable we use 
ECa=1.7EW 8 542. These coefficients are derived empir- 
ically from high S:N spectra where the absolute values 
of [Fe/H] are well known. We test these variations of 
ECa by applying them to our high S:N co-added spectra 
where all three lines are well resolved, such as And XXI 
and XXV, and we find that all three formulae produce 
consistent values of [Fe/H]. The composite spectra for 
each satellite are shown in Figs. [H] and [2U In all cases, 
we find that our results are consistent with photometric 
metallicities derived in previous works. 

In the MW, it has been observed that the average 
metallicities of the dSph population decrea s e wit h de- 
creasing luminosity (e.g.. iKirbv et all 120081 120111) . In 
Fig. [231 we plot the spectroscopic metallicities of the M31 
dSphs (blue dSphs) as a function of absolu te magnitude. 
We p l ot the MW dSphs as red triangles dMartin et all 
20071 IKirbv et all [20081 12011 IBelokurov et al.l 120091: 
Kochet all I2009D . We also include those M31 dSphs 
for which only photometric measure ments of [Fe/H] are 
available (And I, II III, VII and XIV. IKalirai et~aU[20Tol: 
iTollerud et aLll201~2j ). and these are highlighted as encir- 
cled blue points. The dashed line re presents the best-fit 
to the MW dSph population from IKirbv et all (|2011l) . 
The three relatively metal rich ( [Fe/H] ~ -1.5 to -2.0), 
faint (L ~ 1000 L©) MW points are the three ultra faint 
dSphs Willman I, Bootes II and Segue 2 , and t hese three 
were not included in the IKirbv et al.l (|201lD analysis, 
where the best fit MW relation was determined. We 
see that the metallicities for a given luminosity in the 
M31 dSphs also loosely define a relationship of decreasing 
metallicity with decreasing luminosity, and they agree 
with that defined by their MW counterparts within their 
associated uncertainties. However, it is also noteworthy 
that for dSphs with L < 10 6 L Q , the Andromeda satel- 
lites are also consistent with having a constant metallic- 
ity of ~ —1.8. The same levelling off of average metal- 
licity at lower luminosities was noted by IMcConnachid 
(|2012l) . where they note that this break occurs at the 
same luminosity as a break in the luminosity-surface 
brightness relation for faint galaxies. As such, it could 
imply that the denisty of baryons in these systems, rather 
than the total number of baryons, could be the most im- 
portant facor in determining their chemical evolution. 
The error bars we present here are still significant, so 
it is hard to fully interpret this result, but the hint of 
a metallicity floor in these lower luminosity systems is 
intriguing. 

In Fig. [231 we highlight the positions of our three kine- 
matic outliers, And XIX, XXI and XXV, and we see 
that they fall almost exactly on the MW relation. In 
this figure, systems that have experienced extreme tidal 
stripping would move horizontally to the left, as their lu- 
minosity would gradually decrease as stars are stripped, 
but their chemistry would remain unaffected. One would 
expect to see such behaviour only after the stellar com- 
ponent began to be removed in earnest, after the major- 
ity of the dark matter halo had been removed. If their 
central densities were lowered by some a ctive feedback 
mech anism, such as SNII explosions (e.g.. iZolotov et all 
I2012I) . one would expect the objects to become more en- 
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Fig. 21. — Figures of the S:N weighted composite spectra for each dSph. The average metallicity for each dSph, plus the associated 
uncertainty is printed on each panel (continued over page) 
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Fig. 22. — Figures of the S:N weighted composite spectra for each dSph (excluding And XV and And XVI). The average metallicity for 
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Fig. 23. — Spectroscopically d erived [Fe/H] vs. l uminosity for all 
MW (red triangles, taken from Kirby ct al. 201 lj, with additional 
measurements taken from Martin ct al. 2007; Bclokur ov et al] 
[20091 : iKoch et aljl200Sft and M31 dSphs (blue circles, this work). 
The solid line represents the b est fit relati o nship between these 
two parameters as taken from Kirby ct al. (201 J). The dashed 
lines represent the la scatter about this relationship. We see that 
the M31 dSphs follow this relationship very well within their as- 
sociated uncertainties. As we discuss in § \E\ those galaxies with 
L < 10 6 Lq are also consistent with having a constant metallicity, 
which could indicate a metallicity floor in these fainter systems. 



riched and perhaps brighter, moving them up and to the 
right, potentially allowing them to remain on the MW 
relation. To confirm that this was the case for And XIX, 
XXI and XXV, we would require more information on 
the abundances of these objects and their star formation 
histories, which we do not currently possess. 

7. CONCLUSIONS 

Using new and existing spectroscopic data from the 
Keck I LRIS and Keck II DEIMOS spectrographs, we 
have homogeneously derived kinematic properties for 18 
of the 28 known Andromeda dSph galaxies. Using a com- 
bination of their g — i colors, positions on the sky and 
radial velocities, we determine the likelihood of each ob- 
served star belonging to a given dSph, thus filtering out 
MW foreground or M31 halo contaminants. We have 
measured both their systemic velocities and their velocity 
dispersions, with the latter allowing us to constrain the 
mass and densities within their half-light radii. For the 
first time, we confirm that And XVII, XIX, XX, XXIII, 
XXVI and Cass II are dark matter dominated objects, 
with dynamical mass-to-light ratios within the half-light 
radius of [M/L] half > 10 Mq/ Lq. 

For And XXV, a bright M31 dSph (M v = -9.7) 
we measure a mass-to-light ratio of only [M/L] half = 
10.3lg'y M Q / L Q from a sample of 26 stars, meaning that 
it is consistent with a simple stellar system with no ap- 
preciable dark matter component within its la uncer- 
tainties. If this were confirmed with larger datasets, it 
would prove to be a very important object for our un- 
derstanding of the formation and evolution of galaxies. 



We compare our computed velocity dispersions and 
mass estimates with those measured for MW dSphs, and 
find that the majority of the M31 dSphs have very sim- 
ilar mass-size scalings to those of the MW. However, we 
note 3 significant outliers to these scalings, namely And 
XIX, XXI and XXV, who possess significantly lower ve- 
locity dispersions than expected for their size. These 
results buil ds on the ident i ficatio n of three potential out- 
liers in the lTollerud et all (|2012D dataset (And XIV, XV 
and XVI) . We suggest that the lower densities of the dark 
matter halos for these outliers could be an indication that 
they have encountered greater tidal stresses from their 
host over the course of their evolution, decreasing their 
masses. However, these bright systems still fall on the 
luminosity-metallicity relation established for the dSph 
galaxies of the Local Group. If these objects had under- 
gone significant tidal disruption, we would expect them 
to lie above this relation. As such, this remains puzzling, 
and requires dedicated follow up studies to fully map out 
the kinematics of these unusual systems. 

We measure the metallicities of all 18 dSphs from their 
co-added spectra and find that they are consistent with 
the established MW trend of decreasing metallicity with 
decreasing luminosity. 

This work represents a significant step forward in un- 
derstanding the mass profiles of dwarf spheroidal galax- 
ies. Far from residing in dark matter halos with identi- 
cal mass profiles, we show that the halos of these objects 
are complex, and differ from one to the next, with their 
environment and tidal evolution imprinting themselves 
upon the dynamics of their stellar populations. The An- 
dromeda system of dSphs presents us with an opportu- 
nity to better understand these processes, and our future 
work will further illuminate the evolutionary paths taken 
by these smallest of galaxies. 
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APPENDIX 

A. TESTING THE MEMBERSHIP PROBABILITY ALGORITHM 
A.l. The inclusion of a velocity term in the calculation of Pi 

In T12, the authors do not impose a velocity probability criterion for their membership calculations. Instead they 
require all member stars to have a total probability, based on their positions and colors, of Pmember > 0.1, and then 
apply a 3c clipping to this final sample to prevent any outliers from significantly inflating their calculated velocity 
dispersions. In our analysis, we have avoided making any hard cuts to our sample by also utilising prior information on 
the velocities of our expected contaminant populations and member stars. In § 13.31 we tested our velocity probability 
criteria was not overly biasing our final measurements of velocity dispersion to artificially lower values with the 
introduction of an extra parameter, 77, that allows us to add additional weight to stars in the tails of the Gaussian 
velocity distribution. However, we can further test our velocity criterion by removing it entirely from the probabilistic 
determination, and instead implementing the same cuts presented by T12. This involves cutting stars where Pi < 0.1 
(as determined from Pcmd and Pdist), and also by iteratively removing all stars that have velocities that do not lie 
within 3<7 of the mean of the remaining sample. In Table IA"TI we present the results of this on our measured values of 
v r and a v for our full sample of dSphs. For all objects (bar And XXVII, which is a unique case, as described in § 14. lip 
the systemic velocities derived are within ~ 2 — 3kms _1 of one another. The velocity dispersions we measure from 
our full algorithm tend to be slightly higher on average, and this is to be expected as we do not cut any stars from our 
analysis, and therefore outliers in the velocity profile may be assigned non-negligible membership probabilities that 
will allow them to increase this measurement. By and large, these differences are not significant, with the final values 
agreeing to well within their la uncertainties. 

It is interesting that our algorithm appears to perform better when dealing with dSphs where the number of member 
stars is low. This is best demonstrated by And XI (see Fig. IA1|) . Our algorithm identifies 5 stars with non- negligible 
probabilities of membership, clustered around v r ~ — 430 km s -1 . Our full algorithm measures a systemic velocity of 
v r = — 427.5l3'4kms _1 and a velocity dispersion of a v = Lb^'s kms -1 . One of these stars is slightly offset from the 
other 4 with a more negative velocity of v r = —456.8 km s . Although this star has a reasonably high probability of 
being a member based on its distance from the centre of And XI, and its position in the CMD, it does not survive 
the 3tr velocity clipping procedure of T12. As the number of member stars is so low, cutting one star from the sample 
can have a significant effect, and as such, while the T12 procedure determines a very similar systemic velocity of 
v r = —425.0 ± 3.1kms _1 it is unable to resolve a velocity dispersion. This effect is also seen other systems (such as 
And XIII, XVII, XXII and XXVI), although it is typically less pronounced. 

Another regime where our algorithm performs better than that of T12 is where the systemic velocity of the system 
in question is within the regime of the contaminating Milky Way K-dwarfs. An example of this is the unusual system, 
And XIX, where our algorithm measures a systemic velocity of v r = —111. 4 kms^ 1 and a velocity dispersion 
of a v = 4.7ti'4- However, the procedure of T12 is less able to resolve the kinematics of the system, measuring 

v r = —109.3 ± 5.3kms _1 and a velocity dispersion of a v = 1.8~^i'\. The much larger uncertainty on the dispersion is 
a result of including Milky Way contaminants in the sample which can be difficult to cut out without applying prior 
knowledge of the velocity profile of this population. 3a clipping allows outliers to contribute more significantly to the 
measured profile in this instance, increasing the uncertainty. A similar effect is seen in the And XXIV and And XXX 
(Cass II) objects, which also have systemic velocities in the Milky Way contamination regime. 
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Fig. Al. — As Fig. [6] but for And XI. Using our full probabilistic algorithm (including our P ve i term), 5 stars are identified as likely 
members. When implementing the 3<r clipping of T12, the star with v r ~ — 450kms — 1 is removed from the analysis, despite having a high 
probability of membership based on its color (Pqm d) a n d distance from the centre of the dSph (Pdist)- This has a substantial effect of 
the calculated velocity dispersion, as shown in Table |A"T1 



These results lead us to conclude that the inclusion of a P ve i term in our analysis allows us to more effectively 
determine the true kinematics of the systems we are studying. Further, as no cuts to the sample are required using 
this method, it allows for a more unbiased study of the kinematics of dSphs than that of T12. 
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TABLE Al 

Measured velocities (v r ) and dispersions (a v ) for our full sample of 
dSphs as calculated using cuts on probability of membership calculated 

from only color and position of stars, combined with 3(7 clipping on 
velocities vs. those calculated using our full probabilistic algorithm. 





Object 


v & cuts 




Full algorithm 








« r (km s ) 


a v ( km s ) 


» r (km s — x ) 


cr„ ( km s 1 ) 


And 


V 


-391.1 ± 2.9 


10.8+!;° 


-391.5 ± 2.7 


12 - 2 ±i'.9 


And 


VI 


-339.0 ± 3.0 


ii.9ijj;jj 


-339.8 ± 1.8 


- ^ . -4-1 ' R 


And 


XI 


-425.0 ± 3.1 


g+J.a 


-427.5^34 


7-612.8 


And 


XII 


-558.8 ± 3.7 


0+6-8 


-557.1 ± 1.7 


0+4.0 


A r,A 


"VTTT 


90*3 8 -t- 8 A 


g+16.2 


— zu^±.o zjz ^t.y 


n n+8-i 


And 


XVII 


-260.0± 7 ^ 


i-slii 


-254.3 + f 7 


2 9+ 50 
z - a -2.9 


And 


XVIII 


-345.1 ±3.3 


+4.4 


-346.8 ± 2.0 


Q+2.7 


And 


XIX 


-109.3 ±5.3 


i 5 +6.8 
•■•• -1.6 


-111.6±i-2 


4 7+1-6 
4 - '-1.4 


And 


XX 


-454.6±« 


7 7+8.4 
'•'-3.9 


-456.2l3;i 


7 i+3.9 


And 


XXI 




n tj+2.3 


-362.5 ±0.9 


4-51^ 


And 


XXII 


-131.4 ±2.7 


0+3.1 


-129.8 ±2.0 


2.8±i;S 


And 


XXIII 


-236.9 ±2.1 




-237.7 ± 1.2 


7.1 ± 1.0 


And 


XXIV 


-129.2 ±3.6 


0+6.1 


-129.9±J| 


o c+6.6 
a,0 -8.5 


And 


XXV 


— 107.7±i;i 


n + 2.2 

°-°-1.8 


-107.8 ± 1.0 


3.0+1? 


And 


XXVI 


-264.1 ±4.5 


0+4.8 


-261.7±3;| 


8.7+li 


And 


XXVII 


-517.6±gi 


19.3±JS 




14.81" 


And 


XXX (CassII) 




14-llS 9 


-139.8+.^ 





TABLE A2 

THE EFFECT OF INCLUDING LOW S:N DATA IN OUR ANALYSIS ON DETERMINING V r AND <T, 



Object 


AT, 


And XXI 

v r ( km s^ 1 ) 


<r„ ( km s 1 ) 


N* 


And XXIII 

ti^kms" 1 ) 


a v ( km s 1 ) 


AT* 


And XXV 
« r (kms _1 ) 


<7„ ( km s 1 ) 


S:N> 2 
S:N> 3 
S:N> 4 
S:N> 5 


20 
11 

5 
3 


-362.9 ±0.9 
-364.2 ± 0.9 
-363.9 ± 1.5 
-362.6±|J 


o c+0.9 
°-°-0.7 

3.1 ±0.8 
4.0+H 

4 5+2-6 


22 
10 
5 
2 


-238.0 ± 1.2 
-238.3 ± 1.4 
-239.4 ± 1.1 
-239.5 ± 1.9 


6.6 ±1.1 

0.1_ 1 2 
5 7+1.5 
°-'-1.3 

0.0+ 49 


27 
24 
16 
13 


-107.7 ±0.9 
-107.7 ±1.0 
-108.2 ± 1.2 
-109.0 ±1.2 


2.7±1.1 
3.0 ± 1.2 

z.o±\i 

2.8tH 



Note. — xxx 



A. 2. The effect of low signal-to-noise data on measuring v r and o~ v 

o 

For our brightest targets (i < 22.5) the S:N of our spectra is typically > 3 A -1 . However, as our targets become 

o 

fainter, so too their S:N falls. For spectra with S:N> 1.5 A -1 , our pipeline is still able to measure velocities based on 
the Ca II triplet, with reasonable measurement uncertainties. However, it is prudent to check whether the inclusion 
of these velocities, calculated from significantly noisier spectra, has a detrimental effect on our ability to measure the 
kinematic properties of our dSph sample. 

Such a test is straightforward to implement. We have a number of dSphs within our sample for which our probabilistic 
analysis identifies ~ 30 likely members (such as And XXI, XXIII and XXV). We can therefore use these samples to 
impose S:N c uts on our data to see the effect of this on our measurements of v r and a v . We present the results of this 
test in Table rA"2l and our finding is that, as the level of our imposed S:N cut increases (and so the number of included 
stars decreases), the systemic velocity remains more or less constant. The measured velocity dispersion, however, 
shows some variation. In the case of And XXI and XXV, the dispersion increases with increased S:N, however not 
significantly. In both cases the dispersion calculated from the higher S:N data lies well within la of that calculated 
from the lower S:N data. Intuitively, this makes sense as the spectra with higher S:N are likely to have lower velocity 
uncertainties, and so our maximum-likelihood analysis will attribute more of the spread in measured velocities to an 
intrinsic dispersion, rather than to our measurement errors. In the case of And XXIII, we find the opposite to be true. 
As our S:N cut increases, we find that our measured dispersion decreases. This may be because the number of member 
stars in subsequent quality cuts drops off more rapidly for And XXIII than And XXI and XXV. This suggests that we 
should be extra cautious when interpreting our measured velocity dispersions for dSphs where both the average S:N 
of member stars, and the number of member stars, is low. 
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TABLE A3 

The effect of varying sample size (JV*) on the measurements of systemic VELOCITY, V r , AND 

VELOCITY DISPERSION, a v . N STARS FROM AND XXI, XXIII AND XXV ARE RANDOMLY SELECTED, AND THEIR 
PROPERTIES ARE DERIVED USING OUR FULL PROBABILISTIC ANALYSIS. THIS WAS REPEATED 1000 TIMES. VALUES 
REPORTED BELOW ARE THE AVERAGES FROM 1000 REALIZATIONS, WITH THE UNCERTAINTIES REPRESENT THE 
STANDARD DEVIATION OF THE 1000 REALIZATIONS. 



Object 


And XXI 




And XXIII 




And XXV 




JV» 


v r ( km s _1 ) 


cr v ( km s 1 ) 


» r (km s _1 ) 


cr v ( km s 1 ) 


» r (kms _1 ) 


cr v ( km s 1 ) 


4 


-362.6 ±3.1 


5.0 ±3.3 


-237.2 ±3.8 


7.5 ±3.6 


-108.8 ±2.4 


5.4±4.3 


6 


-363.0 ±2.0 


4.0 ± 1.9 


-237.3 ±3.0 


6.9 ±2.6 


-108.0 ±2.7 


4.0 ±2.3 


8 


-362.9 ±2.1 


4.2 ± 2.0 


-236.8 ±2.8 


6.9 ±2.3 


-107.6 ±2.0 


2.7 ±1.1 


10 


-362.7 ± 1.8 


4.1 ± 1.8 


-237.3 ±2.5 


6.5 ±2.1 


-107.9 ± 1.5 


3.1 ± 1.0 


15 


-362.9 ±1.5 


4.3 ± 1.2 


-237.4 ±2.0 


6.8 ±1.4 


-107.9 ± 1.5 


3.1 ± 1.0 


20 


-362.9 ±1.3 


4.3 ± 1.2 


-237.6 ± 1.7 


6.9 ±1.4 


-107.8 ± 1.2 


3.0 ±0.9 


25 


-363.0 ±1.0 


4.4 ± 1.0 


-237.4 ± 1.4 


7.1 ± 1.7 


-107.8 ± 1.0 


3.1 ±0.8 


30 


-362.8 ±1.1 


4.4 ± 1.1 


-237.7 ± 1.3 


7.2 ±1.4 






Full sample 


-362.5 ± 0.9 


*-°— l.i 


-237.7 ± 1.2 


7.1 ± 1.0 


-107.8 ± 1.0 


8.0+" 



A. 3. The effect of small sample sizes on determining kinematic properties of dwarf galaxies 

Obtaining reliable velocities for member stars of faint and distant systems is a difficult task that can only be achieved 
with the largest optical telescopes, such as Keck. Given the demand for facilities such as this, any observing time 
awarded must be used as effectively as possible, and this often means compromising between deep pointings for a few 
objects, and shallower pointings for a number of objects. With longer or multiple exposures on a single target, one can 
build up impressive samples of member stars for an individual dSph. For exa mple, the SPLASH c ollaboration observed 
a total of 95 members in And II, one of the brightest M31 dSph companions (jKalirai et al.ll2010l ) by taking 2 separate 
exposure fields over this large object. However, multiple exposures such as these produce diminishing returns as you 
move down the luminosity scale to fainter, more compact dSphs. This is both because of their smaller size with respect 
to the DEIMOS field of view, and the fewer number of bright stars available on the RGB to target. In this case, the 
only way to identify more members is by integrating for longer, but given the paucity of stars, the trade-off between 
time spent exposing and additional members observed can be quite expensive. Such difficulties inevitably lead to the 
inference of dynamical properties for an entire system from a handful of stars. It is important for us to understand 
the effect this bias has on our results, and how reliable the quoted values are. We test this using our datasets for 
which we identify > 25 member stars, namely And XIX, XXI, XXIII and XXV using the following method. We select 
4, 6, 8, 10, 15, 20, 25 and 30 stars at random from each dataset and then measure the systemic velocity and velocity 
dispersion using our probability algorithm. This was repeated 1000 times for each sample size. In cases where the 
algorithm is unable to resolve a velocity dispersion, we throw out the result and resimulate, as null results here will 
affect our averages and will not inform us whether the instances in which we are able to resolve a veloc ity d ispersion 
from small numbers of stars are producing valid, reliable result. We display the resulting values in Table lA3| with the 
true value recovered from the full sample shown in bold in the final row for comparison. We show that in all these cases, 
the systemic velocity and velocity dispersion are recovered well within the scatter of the 1000 simulations even when 
dealing with sample sizes as small as 4 stars, so long as the measurement is resolved. In cases where we are unable to 
resolve a dispersion, we find that our resulting uncertainties are not meaningful. This is shown explicitly in the case of 
And XVIII, where we can compare our upper limit for the velocity dispersion as determined from our algorithm with 
the dispersion calculated in T12 from a much larger dataset. We see that our uncertainty is not consistent with their 
result. As such, we advise that in all cases where we calculate velocity dispersions from small samples (N* < 8, And 
XI, XII, XX XXIV and XXVI), the dispersion measurements should be treated as indications of the likely dispersion, 
and need to be confirmed with follow-up studies. 

A. 4. Testing our algorithm on the SPLASH sample of M31 dSphs 

In T12, the authors reported on the kinematic properties of 15 M31 dSphs, And I, III, V, VII, IX, X, XI, XII, 
XIII, XIV, XV, XVI, XVIII, XXI and XXII, and the positions, and measured velocities (plus uncertainties) for each 
likely member stars were published as part of that work. The authors were kind enough to also give us access to 
these properties for their non-member stars so that we might run our algorithm over the full samples to see if we 
reproduce their results. Our technique for assigning membership probability differs from theirs in that we use the 
velocities of stars as an additional criterion for membership, whereas they use a cut on both the resulting probability 
(P(member > 0.1), and a 3a clipping on the velocity. In addition, where we use PAndAS CFHT MegaCam g— 
and i — band photometry for our membership analysis, the SPLASH team use their own Washington-DD051 filter 
photometric dataset, in membership classification. As such, sm all d ifferences might be expected, but if our technique 
is robust our results should well mirror those of T12. In Table lA4l we compare our calculated values of v r and a v to 
those published in T12. As the measurements made in T12 for And XI and XII are made from only 2 stars, we do 
not include these in this test. In general, the results from both analyses agree to within la of one another, with the 
majority of them being well within this bound. Typically, we find that our procedure measures slightly larger values 
for a v than that of T12 (with the exception of And IX, And XVIII and And XXII). This is to be expected, as we do 
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TABLE A4 



Object 


T12 analysis 
v r ( km s J 


a v ( km s J 


Our analysis 
» r (kms J 


a v ( km s J 


And 
And 
And 
And 


I 

III 

V 

VII 


-376.3±2.2 
-344.3±1.7 
-oy I .o±l.D 
-307.2±1.3 


10.2±1.9 
9.3±1.4 

lU.Oztl.l 

13. 0± 1.0 


-376.3 ±1.3 
-344.2 ± 1.2 
— oyD.U ± l.U 
-307.1 ± 1.1 


10.6 ± 1.0 
10.1 ± 1.9 

11 A 1 1 O 
11.4 ± l.z 

13.1 ±0.9 


And 


IX 


-209.4±2.5 


10.9±2.0 


-210.3 ± 1.9 


10.2+J"? 


And 


X 


-164.1±1.7 


6.4±1.4 


-165.3 ±1.5 


6-o±i-.i 


And 


XIII 


-185.4±2.4 


5.8±2.0 


-183.0l£* 


S.6+?* 


And 


XV 


-323±1.4 


4.0±1.4 


-322.6 ± 1.1 


e.oi r ;° 


And 


XVI 


-367.3±2.8 


3.8±2.9 


-366.il*;? 


4 2+ 4 - 8 


And 


XVIII 


-332.1±2.7 


9.7±2.3 


-330. 7±|-_? 


r.5l|;f 


And 


XXI 


-361.4 ±5.8 


7.2 ± 5.5 


-358.91^ 


8.5l 6 5 ;? 


And 


XXII 


-126.8 ±3.1 


o 54+4.16 


-124.2l| : | 


o.o+ 5 - 7 



not cut stars from our analysis based on their velocity, instead we down-weight their probability of membership. As 
such, those stars considered as outliers would naturally inflate our dispersions above those measured by T12, but the 
effect is marginal. These results demonstrate that our technique for assigning probability of membership of individual 
stars within M31 dSphs based on their photometric properties and velocities is robust, and comparable to that of T12. 
However, as discussed in § IA.11 we find that our technique is superior as it requires no cuts to the final dataset to be 
made, reducing the bias in these measurements. 

B. REANALYZING OUR PREVIOUSLY PUBLISHED RESULTS 

To ensure our analysis of the global properties of the M31 dSph population in § 15.11 is homogenous, we reanalyzed 
our previously published datascts using our new probability algorithm. Here we briefly summarize the results of this 
analysis and compare the new results to the published works. The objects we discuss here are And V, VI (first 
published in lCollins et"aI|[20TU) . XI, XII and XIII (published in lChapman"et~^[2005l [20071: iCollins et al.ll2010D . And 
V was observed with the LRIS instrument, while the remaining objects were observed with the DEIMOS instrument. 

B.l. Andromeda V 

Andromeda V (And V) was observed using the LRIS instrument on Keck I rather than the DEIMOS instrument on 
Keck II. LRIS has a lower resolution than DEIMOS, and a smaller field of view, which lowers the accuracy of velocity 
measurement and limits us to only ~ 50 targets within a mask compared with 100 — 200 for a DEIMOS mask. The 
raw data were also not reduced using our standard pipeline, owing to problems with the arc-lamp calibrations, and 
were instead analyzed using the NOAO.ONEDSPEC and NOAO.TWODSPEC packages in IRAF. The results from 
this reduction, plus an analysis of the data using hard cuts in velocity, distance and color to determine membership 
were first published in Collins et al.l (|2011h . Our full probabilistic analysis identifies 17 stars with a non- negligible 
probability of belonging to And V. Our technique determines a systemic velocity of v r — —391.5 ± 2.7 km s -1 and 
a v = 12.2l?'g kms -1 . Comparing these values to our previously published results (v r — —393.1 ± 4.2kms~ 1 and 
a v = ll.Sl^'gkms -1 , ICollins et alJ[20Tll ) we find them to be consistent within the quoted uncertainties. We also 
compare our results to those of T12, who measured v r = 397.3 ± 1.5 km s -1 and a v — 10.5 ± l.lkms -1 from a larger 
sample of stars (85 members cf. 17) using the higher resolution DEIMOS spectrograph. The velocity dispersions of 
both are consistent within their la uncertainties, as are the systemic velocities. Given the difference of a factor of 5 in 
number of probable member stars between our study and that of T12, this consistency is reassuring, and demonstrates 
the ability of our technique to accurately determine the kinematics of M31 dSph galaxies from small sample sizes. 

B.2. Andromeda VI 

As And VI sits at a large projected distance from the centre of M31 (~ 270 kpc), it was not observed as part of 
the PAndAS survey and we were una ble to use CFHT d ata for our Pcmd determination. Instead, we use Subaru 
Suprime-cam data (PL N. Arimoto, see lCollins et al.ll2011l for a full discussion of this data). Using our full probabilistic 
analysis, we identify 45 stars with Pi > 10 -6 . Our technique determines a most-likely v r = —339.8 ± 1.8 km s -1 , 
with a v = 12.41^3 kins" 1 . Comparing these values with the results of ICollins et al.l (|2011| ) who measured v r = 

-344.8 ±2.5 km s" 1 and a v = 9.4!^ kms" 1 , we see that both the systemic velocities and the velocity dispersions are 
consistent within quoted uncertainties. The slight differences between our previous study and this work are simply a 
result of the application of our new technique. 

B.3. Andromeda XI 

The kinem a tic pr operties for Andromeda XI as measured from this DEIMOS data set were first published in 
ICollins et al.l (|2010l ). Here we identify 5 stars as probable members. We determine most-likely parameters of 
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v r = -427.5^3:4 km s" 1 and a v = 7.6±%% kms -1 . In iCollins et al.l (|2010T > we measured v r = -419.4j:f;|kms _1 
and we were unable to resolve the velocity dispersion for the dSph, measuring a v = O.O+^kms" 1 (where the upper 
bound represents the formal la uncertainty on the unresolved dispersion), which implied a higher systemic velocity 
and lower velocity dispersion for And XI. However, in that analysis, one star with a velocity of ~ — 440 kms -1 was 
considered to be an outlier based on its velocity, and thus excluded from the kinematic analysis. Here, our algorithm 
gives this star a non-zero probability of membership, which likely decreases the systemic velocity and increases the 
dispersion. T12 also presented observations for And XI, but they were not able to cleanly detect the galaxy. They 
identified 2 stars with highly negative velocities (~ — 460 kms -1 ), which are offset from our systemic velocity by 
~ 30 kms -1 . Given the very negative velocities of their 2 stars, the probability of them both being M31 contaminants 
seems low, and some other explanation may be more suitable. Between observations, there is one star common to both 
(a =00:46:19.10,(5 =+33:48:4.1), for which we measure a velocity of —427.16 kms -1 compared with —461.6 kms -1 . 
This amounts to a statistical difference at the level of 5a. One obvious avenue to check is that there has been no 
velocity offset introduced by a rogue skyline that falls within the region of one of the three Ca II lines. We have 
carefully checked the spectra of each of our 5 probable members to see if this has been the case. We also rederive 
the velocity based on cross-correlations with each of the three lines individually, rather than with the full triplet. We 
find these results, and their average to be entirely consistent within the associated errors from the velocities derived 
using the technique discussed in § 12.21 This large discrepancy is puzzling, particularly as the methods used to measure 
velocities in this work are almost identical to those of T12. Without further data it is not possible for either team 
to pin down the exact issue, or which of the datasets gives the true systemic velocity. This argues for taking further 
observations within these faint systems, in the hopes of better understanding both the systems themselves, and any 
systematics introduced by the DEIMOS instrument. It is comforting that, in all cases where this offset is observed, 
the velocity dispersions measured by each team are consistent with one another, suggesting the problem affects all 
observations identically. 

In this case, we have identified a greater number of potential members stars than in the work of T12, the spectra for 

o 

all of which have relatively high S:N (S:N^ 5 — 10 A -1 ). Therefore, our results should be considered as more robust 
than those of T12. 

B.4. Andromeda XII 

Kinematic properties for Androm eda XII as determined fr om the DEIMOS data set presented here were previously 
published in iChapman et all ()2007f ) and ICollins et al.l (|2010l ). and in both cases, membership was largely determined 
using hard cuts in velocity. This object has been of particular interest as it possesses an extrem ely negative systemic 
velocity with respect to Andromeda, suggesting that it is on its first infall into the local group fsee lChapman et al.ll2007l 
for a full discussion). Using our new algorithm, we measure v r — —557.1 kms -1 and an unresolved velocity dispersion 
of o~ v = 0.0 +40 kms -1 , where the upper bound represents the formal la uncertainty on the measurement. As the 
dispersion is unresolved, the lower error bound is undefined , as having a ne g ative veloci ty dispersion i s unp hysical. 
These values are completely consistent with the results of IChapman et al.l (|2007l ) and ICollins et al.l (|2010l ) (v r — 
— 558.4 ± 3.2 kms -1 and a v = 2.6^26 kms -1 ). T12 also presented observations for And XII, but as for And XI, they 
were not able to cleanly detect the galaxy. They identified 2 stars with highly negative velocities (~ —530 kms -1 ), 
which are offset from our systemic velocity by ~ 30 kms -1 . In this instance, both the stars they observed overlap with 
two of our likely members, situated at a=00:47:28.63,(5=+34:22:43.1 and a =00:47:24.69,(5 =+34:22:23.9, and these 
velocities are offset from those that we measure at a statistical level of 3.8a. This suggests that the self-same calibration 
effect that causes an offset between our results for And XI and those of T12, is present here also. In the case of And 
X II, our mask was obser ved over two separate n ights, giving us two velocity measurements for each star (as discussed 
in IChapman et al.l 120071 and ICollins et al.l 120101 ) , and we saw no evidence for systematic offsets of this magnitude in 
the night to night velocities, making a calibration error within our dataset seem unlikely, though not impossible. We 
therefore conclude that, owing to our larger sample of members and repeat observations, our measurements for the 
kinematic properties of And XII are more robust than those of T12. 

B.5. Andromeda XIII 

The kinematic properties of Andromeda XIII were also presented in ICollins et al.l ()2010[ ). And XIII sits at a large 
projected distance from Andromeda ( ~ 120 kpc) in the southern M31 halo, so we expect contamination from 
the Milky Way and Andromeda halo to be low. It is surprising then, that we see significant structure within our 
DEIMOS field. This is also seen within the 3 fields observed by T12, who attribute this over-density of stars located 
at v r ~ —120 kms -1 to an association with the TriAnd over-density within the Galactic halo. We too see a number 
of stars between —140 kms -1 and —100 kms -1 . From their positions within the CMD of And XIII, they appear more 
consistent with MW foreground K-dwarfs than M31 RGB stars. As such, these are also likely associated to this MW 
substructure. 

Using our full probabilistic analysis, we identify the most probable And XIII stars as those 4 that cluster around 
v ~ —200 kms -1 , and we determine v r — —204.8 ± 4.9 kms -1 and are unable to resolve a velocity dispersion, with 
a v = 0.0 +81 , where the upper limit indicates the formal la uncertainty on the likelihood distribution. Given the very 
large uncertainties on these values (mostly a factor of the low number of member stars) it comes as no surprise perhaps 
that these results are consistent with the results in ICollins et al.l (|2010D (v r = — 195.0tg4kms -1 and a v = 9. Tt|J), 
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although not with those of T12 (v r = — 185.4 ± 2.4 km s -1 and a v — 5.6 ± 2.0 km s" 1 ), where they derived parameters 
from three times the number of member stars that we present here. Given the highly substructured nature of the And 
XIII field, and the fact that our detection is at a very low significance (only 4 stars), we find that the T12 measurements 
for the kinematics of And XIII are more robust than ours. 



